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PREFACE 


The experiment and measurements described in this report were 
performed under the Radio -Frequency Interference and Propagation 
Program (RIPP) administered by the National Aeronautics and Space 
Administration (NASA). The planning for this program was guided by 
a steering committee composed of members (and their technical 
advisors) representing the Federal Communications Commission (FCC), 
the Office of Telecommunication Management (OTM - now the Office of 
Telecommunication Policy), NASA, and the Institute for Telecommuni- 
cation Sciences (ITS). The technical investigations were carried out 
by the ITS in two basic phases that spanned one year from September, 

1970 to October 1971. This report covers only Phase I of the program, 
at which time precipitation scatter data were obtained on frequencies 
in both the S-band and X-band portions of the electromagnetic spectrum. 

The experimental period for Phase I was from late September, 1970, 
through January, 1971, and this report contains the computer analyses 
of data for this four -month period. 

Phase II of the program included a reconfiguration of the experiment, 
that eliminated some of the original transmission paths and the X-band 
measurements. Computer analyses as noted in this report were not continued* 
and a manual analysis was performed by NASA for only the S-band data 
for the balance of the year. Results of this work are reported in a NASA 
document (1972) listed in the references. 

Magnetic tape data for Phase II of the experiment are on file at the 
Institute for Telecommunication Sciences (ITS) laboratories in Boulder, 
Colorado, and are available for additional analyses. 
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AN EXPERIMENTAL STUDY OF THE TEMPORAL STATISTICS OF 
RADIO SIGNALS SCATTERED BY RAIN 

R. W. Hubbard, ,J. A. Hull, P. L,. Rice, and P. I. Wells * 

ABSTRACT 


A fixed-beam bistatic CW experiment designed 
to measure the temporal statistics of the volume reflec- 
tivity produced by hydrometeors at several selected alti- 
tudes, scattering angles, and at two frequencies (3. 6 and 
7.8 GHz) is described. Surface rain gauge data, local 
meteorological data, surveillance S-band radar, and great- 
circle path propagation measurements were also made to 
describe the general weather and propagation conditions 
and to distinguish precipitation scatter signals from those 
caused by ducting and other nonhydrometeor scatter 
mechanisms. The operating characteristics of the various 
system components used in the experiment are presented. 
The data analysis procedures were designed to provide an 
assessment of a one -year sample of data with a time resol- 
ution of one minute. The results to date cover the time 
interval of September 15, 1970 to January 31, 1971. The 
cumulative distributions of the bistatic signals for all of the 
rainy minutes during this period are presented for the sev- 
eral path geometries. These cumulative distributions of 
amplitude -time durations of hydrometeor scatter signals 
are compared with measured cumulative distributions of 
surface rain-rate measurements obtained for the same 
time periods. The surface rain-rate statistics are com- 
pared with long-term excessive precipitation data available 
for the southeastern United States. 

Key Words: Propagation, Precipitation scatter. Rain-rate 
statistics. Coordination distances. Satellite/ 
terrestrial communications. Frequency sharing. 


*The authors are with the Institute for Telecommunication Sciences, 
Office of Telecommunications, U. S. Department of Commerce, Boulder, 
Colorado 80302. 
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1. INTRODUCTION 


As an integral part of the Radio -frequency interference Propa- 
tion Program (RIPP) administered by the National Aeronautics and 
Space Administration (NASA), the Institute for Telecommunication 
Sciences (ITS) conducted an experiment at Langley, Virginia, to 
study the propagation/interference problem caused by scatter from 
precipitation. 

The objective of this fixed-beam, bistatic CW experiment was to 
measure the temporal statistics of the volume reflectivity produced by 
hydrometeors at several selected altitudes and scattering angles, and to 
determine the relationship of these measured path-loss statistics to 
those of local surface rain-rate measurements. Extensive meteoro- 
logical, surveillance radar, and great-circle path measurements were 
also made to describe general weather and propagation conditions. 

When the earth terminals of a satellite communication system 
share frequency bands with terrestrial microwave systems, the ad- 
dition of a new terminal to„either system requires that it be established 
in such a manner that no harmful interference will result to either the 
new terminal or any of the terminals of the existing systems. It is 
necessary, therefore, in the planning, design, and construction of such 
systems to have estimates of expected interference to or from the 
proposed new system. 

Current procedures for the estimation of interference between 
terrestrial radio systems and the earth stations of satellite communi- 
cation systems are based upon reports of the CCIR, especially CCIR 



Report 382 [1966a] as revised in New Delhi. The details of these pro- 
cedures are given in reports by the Federal Communications Commission 
[1965] * and the COMSAT Corporation (unpublished report 19 69). 

Many factors in the design and location of a telecommunication 
system can be optimized to provide a maximum ratio of wanted to un- 
wanted signals. Propagation effects, however, are essentially beyond 
the control of the system designer. They enter into the system per- 
formance by introducing distortion of the desired signal. This may be 
independent of the ratio between the desired and undesired signal [CCIR, 
1966b] (for example, phase interference due to multipath, pulse distortion, 
intermodulation, etc. ) or may involve power fading and interference effects 
that degrade the quality of service by reducing the desired to undesired 
signal ratio. In the latter case, a "safety margin" can normally be applied 
to assure an adequate power level for the victim service. The goal of 
interference probability modeling is to develop a method for determin- 
ing the required "safety margin" as a function of the probability of harmful 
interference that can be determined or predicted. The development of 
models for this purpose has been a contribution of the Institute for 
Telecommunication Sciences (ITS) and its predecessors over the past 
two decades [Rice et al. , 1967]. The methodology and mathematical 
development of such models is based on the concept of service probability 
[CCIR, 1966c], and has been treated in a summer seminar series pre- 
sented by ITS in 1970. 

Since most telecommunications systems can tolerate inter- 
ruption for relatively short, infrequent time periods and must operate 
on essentially a continuous basis in the absence of interference, the role 
of propagation is quite different in performance predictions from its 
role in interference probability predictions. For performance prediction 
models (in the absence of interference) only one dominant propagation 
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mechanism normally needs to be considered. The objective of such a 
model is to predict the probability that this mechanism will be present 
a high percentage of time and be sufficiently efficient to assure that an 
adequate signal will be propagated along the desired path. In contrast, 
the interference probability model that is of concern in this report involves 
several possible propagation mechanisms that can couple interfering signals 
between services sharing the same frequencies. Some of these mechanisms 
may be associated with more than one path. They may be highly efficient 
and may exist for long periods of time. In the performance prediction 
models, the statistical properties of the propagation mechanism are usually 
stated in terms of the distribution of hourly median values, while in the 
case of an interference probability model the statistical properties of rare 
events are needed. 

In considering the aspects of sharing between space and terrestrial 
systems, CCIR Study Group 5 [CCIR, 1970] has identified seven propaga- 
tion factors which should be considered in the development of methods to 
predict the interference between the services, namely: effects of terrain, 
effects of vegetation and buildings, the effects of random tropospheric 
inhomogeneities, the effect of tropospheric ducts or layers, scatter from 
hydrometeors, scatter from aircraft, terrain diffraction and tropospheric 
refraction. Among these factors, one that has been least understood 
in terms of the statistical parameters required to predict the 
interference potential, is scatter from hydrometeors. This particular 
factor is the subject of this study. 

It is not feasible to measure the precipitation reflectivity directly 
in all system configurations that must be coordinated to prevent detri- 
mental rain-scatter interference. Thus, there is a need for a model 
based upon an available data source, such as a meteorological parameter 
that is measured generally on a worldwide basis. The readily available 
data are primarily surface rain measurements. 
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The modeling required for this form of propagation interference 
consists of the selection of pertinent information regarding the temporal 
variations of geophysical parameters, and the electromagnetic wave 
interaction with these parameters, to provide an anticipated statistical 
performance of a given system. A primary concern in this work is 
a precipitation scatter model that will provide an appropriate statisti- 
cal basis for the propagation factors associated with the scatter from 
hydrometeors. The physics of scattering by hydrometeors and turbulent 
layers has been treated and experimentally validated as described in a 
recent report by Crane [l970a] . His work has developed models that 
compare radar backscatter and bistatic scatter phenomena for both 
turbulent layer scattering and scatter by precipitation. These models 
are based upon the reflectivity q associated with the scatter region, 
which was measured directly in the validating experiment with a high- 
power precision radar system. 

The statistics of surface rainfall rates discussed in section 5.3 
and in Rice and Holmberg [1972] are based on close examination of ten 
years of surface data for 49 measuring stations in the U.S. Recorded 
rainfall rates are ratios of the depth of water collected in standard rain 
gauge instruments and the period of time of the record* ranging from 5 to 
180 minutes, and are expressed in mm/hr. In order to estimate the rain- 
fall rates within a shorter period of time (i.e., one minute), extra- 
polation techniques were applied, and a model for the cumulative time 
statistics of one-minute rainfall rates was developed. Two basic "modes" 
are defined* with different average rainfall rates, R and R . The 

X £ 

data base for the model is large enough to describe an "average year" 
(from 1951 - I960) for each of several climatic zones in the U.S. 
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The fixed-beam experiment at Langley, Va. was designed 
to obtain measured rainfall rates with short resolution times, in order 
to directly develop the type of cumulative statistics used in the Rice and 
Holmberg model over the span of a one -year measurement interval. 

In addition, bistatic CW measurements were made of the volume reflectivity 
caused by hydrometeor scattering at several scattering angles. The 
surface rain-rate measurements obtained are reordered into a cumu- 
lative distribution and compared with the statistical prediction for 
an average year, using the Rice and Holmberg model. A long data 
sample will permit an empirical relation between the surface rain-rate 
statistics and the altitude -dependent volume reflectivity statistics to 
be made. Using these relationships, we may then apply the model to 
other locations (climatic regions) to predict the interference probability 
based on surface rain measurements for that location. 

This report discusses the physics of the experiment, some 
statistical limits of rain scatter, and the expected tropospheric path 
loss for the geometric configuration of the experiment. The geometric 
configuration is described, and the system operating characteristics 
required for the evaluation of the observed signals and events are 
presented. The data processing and reduction procedures for both 
the rain gauge and radio path loss measurements are described, 
and the cumulative data for a four -month period of the experiment 
are presented. Finally, some comparative analyses of surface rain- 
rate and radio data are presented to demonstrate the performance of 
the system, and to indicate the level of confidence that can be placed 
in the resultant data. 
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2. DESCRIPTION OF THE EXPERIMENT 
2.1 Physics of the Experiment 

Precipitation scattering of radio energy in the frequency bands 
shared by terrestrial microwave and satellite communications systems 
has been clearly established [Carey and Kalagian, 1969; Gusler and 
Hogg, 1970] as a dominant source of potential mutual interference in 
these systems, for time periods exceeding the present allowable limits 
established by the Federal Communications Commission (FCC). 

The Virginia fixed-beam experiment was designed to provide 
measurements of the path loss variations caused by hydrometeor 
scattering at four discrete altitudes and at two frequencies representa- 
tive of the low and high ends of the present bands allocated for satellite - 
terrestrial communication systems. Surface rain-rate measurements 
were made at the sub-points of the common volumes defined by the inter- 
section of several transmitter antenna beams and the receiver main 
beam (see figure 1, page 28). From these measurements, the long- 
term statistics of the path loss (or volume reflectivity) are accumu- 
lated, as are the corresponding surface rain-rate statistics. Unwanted 
signals can occur in a given communication system configuration by 
means other than hydrometeor scatter; for example, atmospheric 
ducting. Thus, it wag necessary to obtain sufficient ancillary meteoro- 
logical data during this experiment to permit separation of the non- 
rain scatter mechanisms that are possible, assuring that the statistics 
of the path loss measurements are properly interpretable. The design 
of this fixed-beam propagation experiment has emphasized measurements 
that permit proper identification of scatter signals due to rain and hail 
in the common volumes. 
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Most of the physics involved in the experiment, whether meteo- 
rological or electromagnetic in nature, are visible in the bistatic radar 
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is power received in the direction (0,<p) in watts, 
is power transmitted in the direction (0 ',<p 1 ) in 
watts 

are respectively the azimuth and elevation angles in a 
polar coordinate system with their origin at the receiver, 
are the appropriate gain factors for the receiving 
and transmitting antennas relative to isotropic, 
are the ray ranges from the transmitter and 
receiver, respectively, to the element of common 
volume dv in meters, 

is the ray distance from the transmitter to the near 
"edge" of the scattering substance in meters, 
is the ray distance from the receiver to the near 
"edge" of the scattering substance in meters, 
is clear atmosphere absorption (or attenuation) 
coefficient of a beam, 

is the attenuation caused by precipitation and/or 
turbulence, 

is the volume of the scattering region in cubic meters. 
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V 


is the local reflectivity of material in the common 
volume, V, in inverse meters, 
is the wavelength of incident and scattered power 
in meters, and 

denotes the volume of integration. 

The important physical aspect of the experiment is the reflec- 
tivity of energy by hydrometeors in the intersection of the transmitter 
and receiver antenna beams, which determine the common volume, V. 

Thus, the reflectivity factor, T], is the key unknown in equation (1). 
The experiment will distinguish between clear air turbulence and rain- 
fall, although there remains a problem of discrimination in a convective 
precipitation situation where both turbulence and heavy rainfall can co- 
exist with near -forward scattering. The precise reflectivity encountered 
in any given common volume intersection is generally not known. In 
the case of precipitation, the reflectivity is dependent upon the size and 
distribution of the hydrometeors within the common volume. To deter- 
mine 77 precisely requires that measurements be made directly in the 
common volume, for example, the backscatter measurements made with 
a precision monostatic radar in the work noted previously by Crane [l970a]. 
The attenuation and absorption terms in equation (1) may gen- 
erally be neglected at the lower (S-band) frequencies used in this experi- 
ment without appreciable error. At the X-band frequencies, attenuation 
due to rain can be significant if the rain fills a major part of the beam of 
both the receiving and transmitting antennas. However, if the rain is con- 
centrated in and near the common volume (of primary interest here), the 
attenuation may also be neglected at these higher frequencies. 


X 
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The effect of the terrain and obstacles such as buildings or trees 
in the immediate foreground of the antenna can be a significant factor. 
Terrain or other obstacles should be considered when determining the 
effective power radiated from, or received by, an antenna at any specific 
direction [Q Also, the gain factors g^ and g include the effects of 

directivity and polarization of the transmitting and receiving antennas in 
the interference problem. Considerable effort (private communication 
from R. R. Bergman and L. G. Hause, 1970, of the ITS staff) was 
expended toward the determination of the antenna patterns and hence the 
description of the antenna gains in a particular direction for use in (1). 
Receiving antenna pattern measurements were made using a small target 
transmitter located on a tower at Poquoson, Virginia. Representative 
transmitter antenna patterns were measured on the antenna range 
maintained by NASA at Langley, Virginia. Assuming that there is no 
attenuation or angular dependence of the scattered power in the reflectivity 
factor Tj , a simplified version of the bistatic radar equation (1 ) may be 
used, where the antenna gain factors are as determined above. The 
simplified equation is 


P t x s r g t n v 

3 2 2 

64tt sf S 

t r 
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where p and p are now, respectively, the mean power received and 
r t 

transmitted along the main axes of the antenna beams, and S and S 

t r 

are now the range values measured to the intersection of the axes within 
the common volume intersections. The average reflectivity factor V 
may be related to the radar reflectivity value Z by 
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where |k|^ is a parameter dependent upon the dielectric properties 

of the scatter material. In the case of rain and for frequencies between 

2 

1 and 35 GHz, |k| % 1 [Crane, 197 0a], and is taken as 0.93 (an 

approximate value for water). For purposes of this experiment, p 
may thus be approximated by 


- 5-4 

-18 

(4) 

p = 0. 93 TT x 

Z x 10 

6 , 3 

where Z is in the units of mm /m as 

- . . -1 

derived from rain 

distribution 


theory, and p is in m 

For most of the common volumes in the Viginia experiment, 
the volume of intersection between antennas may be approximated by a 
thin cylinder, the intersection of a pencil beam (receiver) with the 
broader beam of the transmitting antennas. In these cases, the volume 
is given approximately by 

V = rr/4 (a S ) 2 a S esc 0 (5) 

r r t t 

where 

a is the half-power beamwidth of the receiving antenna 

in radians, 

a is the half-power beamwidth of the transmitting antenna 

in radians, and 

e is the scattering angle in degrees (the supplement of the 

angle between antenna axes at the scatter intersection. ) 
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Substituting (4) and (5) into (2) yields 


(p, X 2 g g)(0. 93 7T 5 X 4 Z x 10 18 )( 7T /4)(o; 2 a, CSC 0) 
t r t r t 


64 rr 3 S, 


( 6 ) 


Rearranging terms of (6) and adding a factor to account for all line 
losses we may express the total power at the receiver input terminal 
in dB as 

P = P + G + G - L - 20 log X + 1 0 log Z - 1 0 log S 

T t T t t t 

+ 10 log A + 10 log (0.93) + 10 log ( 7T 5 x 10' 18 ) 

- 10 log (64 7 7 3 ) (7) 


where 


10 log p dBm, 

10 log p^ dBm, 

10, log g dB, 
r 

10 log g t dB, 

total line losses in dB, and 
2 

( n / A) a a esc 0, 
r t 


Evaluating terms in (7), and inserting the transmitter power used in the 
experiment, 10 watts (40 dBm), results in 

P = -148.4 + G + G - L - 20 log X + 10 log Z 
r r t t 


- 10 log S^_ + 10 log A. 


( 8 ) 



Two nominal transmission frequencies were used in the Virginia 
experiment as noted earlier. Each of these was subdivided into a number 
of transmitting frequencies as outlined in section 2.3. For the purpose of 
equation (8), the wavelengths may be determined for the nominal frequen- 
cies without significant error in later calculations. Thus the nominal 
frequencies and free space wavelengths are 

f = 3. 672 GHz, 

s 

X = 300/3672 = 0,0817 meters, 

s 

f = 7.834 GHz, and 

x 

X = 300/7834 = 0. 0383 meters, 

x 

The receiving antenna at the Langley, Virginia, earth station (ES) 
was a 30-ft parabolic dish, common to both the S- and X-band frequencies. 
The measured gains of this common-feed antenna were, respectively, 

(G ) = 47. 5 dB at 3. 672 GHz, and 

r s . 

(G ) = 50.8 dB at 7.834 GHz. 

r x 

The transmitting antennas used were 10-ft parabolic dishes at S-band 
and 6-ft parabolic dishes at X-band. The gains of the transmitting 
antennas were, respectively, 

(GJ = 38. 8 dB at 3. 672 GHz, and 
t s 

(G ) = 41.0 dB at 7.834 GHz. 

t x 
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Substituting the above values into (8) yields 


(P ) = -40.74 - L + 10 log Z - 10 log S + 10 log A and (9) 

r s t t 

(P ) = -28. 41 - L + 10 log Z - 10 log S + 10 log A (10) 

r x t t 


for the S- and X-band frequency cases, respectively. The remaining 
nonnumerical terms in (9) and (10) are functions of particular systems 
and path geometry, with the exception of Z. The accepted relationship* 
between the radar reflectivity Z and the rain rate R in the scatter 
volume is given by 


Z = 200 R 


1 . 6 


(ID 


6 3 

where R is in mm/hr, and Z is in mm /m . 
v 

Equation (11) can be converted to logarithmic form as 


10 log Z = 23 + 16 log R v . (12) 

When (12) is substituted into (9) and (10), these equations become 

(P ) = -17.74 - L - 10 log S + 10 log A + 16 log R and (13) 

r s t t v 

(P ) = -5.41 - L - 10 log + 10 log A + 16 log R . (14) 

r x t t v 


#Note: The relation Z = 200 R * appears to be the best empirical 
fit for lower rain-rates. For rain-rates exceeding 75 mm/hr, the 
expression Z = 400 R 1, becomes a better fit. [Austin, 1969] 


13 



(15) 

(16) 


These equations may be further simplified to the form 

( P ) = K +16 log R dBm and 

r s s v 

(P ) = K + 16 log R dBm 
r x x v 

where K and K are constants in dB determined from the individual 
s x 

path parameters. 

The values of these constants have been computed and are given 

in table 1 for all of the bistatic path configurations in the Langley 

experiment. Their values may be viewed as the level of detectable rain 

scatter signal that will be received over each path, if R is assumed to 

be 1 mm/ hr. The individual path and site parameters involved in the 

computations may be found in table 6 of section 3. 

The above development is based on an approximation of the 

common volumes, assuming the intersection of two cylinders as in (5). 

A more accurate computation of the common volumes has been carried 

out by M. J. Miles of the ITS staff (private communication, 1970) using 

conical beam patterns defined by the half-power points of the transmitter 

and receiver antenna beams. The computer-derived values from this 

procedure are also given in table 6, section 3. The difference between the 

two estimates of the volumes is largest for the path configurations associated 

with the Quantico site (D and D , see table 1.). A quantitative comparison 

11 16 

of these values results in less than 1 dB difference in the calculated 
value of the K factors for the Quantico paths. Thus it is concluded that 
the approximation of the volumes as computed from the cylindrical geometry 
of (5) is adequate. 

For the purposes of this experiment, equations (15) and (16) 
may be used in either of two ways to evaluate the measurements. 
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First, the measured interference power at the receiver for any bistatic 

path can be used directly to compute an effective rain-rate (R ) at the 

e 

common volume associated with that path. Second, the measured rain-rate 

at the surface of the earth (R ) may be used to compute the expected scatter 

power if the surface rate corresponds to the actual rate in the common 

volume. Both of these techniques are used in the analysis of data (section 5), 

but not in real or chronological time. For example, it is not anticipated 

that the calculated values of R will correlate closely in real time with the 

e 

measured values of R , as the time correlation between these parameters 
with respect to any particular surface rain gauge instrument will generally 
be destroyed by the non-homogeneous nature of storms, and by wind shears 
which occur during a storm. It is expected, however, that if the values of 
R^ and received power P reordered from chronological time to a 
rank ordering of magnitude values over a reasonably long period, agree- 
ment will be found between the statistics of the two measured values. 

In addition, there must be a height dependence associated with an 
interchangeability between surface rain-rate measurements and effective 
rain-rates or reflectivity measured from common volumes aloft. For 
example, if the antenna beams intersect 30, 000 ft above the surface, very 
little rain in the scatter volume is expected [CCIR 1971]. The fixed- 
beam experiment described here is not designed to determine this height 
dependence in any precise sense. To do so would require a tremendous 
amount of additional data to be recorded directly from the common volume, 
as in the case of the precision radar used by Crane [ 1970a]. The results 
of the fixed-beam experiment will indicate however, some empirical height 
dependence for the particular region of the experiment, for common volumes 
up to 30, 000 ft. It is expected that height dependence models will be de- 
veloped from future experimental and analytical work that is being carried 
out by Crane [1970a], Austin [1971], and Altman [ 1970 ]; one promising 
model has been proposed recently by Dutton and Dougherty [1973], 
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2. 2 Sensitivity of the Experiment 

The maximum sensitivity of the experiment for the parameters 

of interest can be estimated from equations (15) and (16), using the 

values of K and K in table 1 and known characteristics of the 
s x 

receiver. For example, the minimum detectable signals at the 
receiver were approximately 

P (min) = -130 dBm 
r s 

P (min) = -125 dBm 
r x 

for all of the S- and X-band channels, respectively. These values 

can be used as the minimum received power level in equations (15) 

and/or (16), and thus provide an estimate of the minimum detectable 

rain rate. Solving either (15) or (16) for R^, we find, for the minimum 

effective rain-rate R (min), 

e 

_ , . . , -1 r(min) 

R e (min) = log . (17) 

In addition, equations (4) and (11) may be used to compute the 
minimum detectable values of r\ and Z, respectively. The values of 
these parameters for the range of antenna/path configurations of table 1 
are presented in table 2. Note that the values are based upon the 
assumption that each common volume intersection is completely filled 
by precipitation. 


17 



Table 2. Sensitivity Values of the Experiment 


System 
& Path 

R (min) 
e 

mm/hr 

Z (min) 

6, 3 
mm / m 

7? (min) 

-1 

m 

S-band 







- 10 

°11 

0. 47 

58. 8 

3. 76 x 10 

E n 

0.47 

58. 8 

3. 76 x 10" 10 

c 22 

0.08 

3.6 

2. 32 x 10 -11 

C 33 

0.43 

51. 3 

3. 27 x 10“ 10 

D 33 

0.43 

52. 5 

3. 35 x 10“ 10 

B 45 

0. 37 

39.8 

2. 54 x 10“ 10 

C 45 

0. 37 

39.8 

2. 54 x 10" 10 

D 44 

0. 55 

75.8 

4.84 x 10“ 10 

E 4 5 

0. 37 

40. 7 

2. 60 x 10“ 10 


X-band 


D li 

0.40 

45. 7 

6. 03 x 10“ 9 

E 16 

0.39 

43. 6 

5. 75 x 10“ V 

Q 

C 37 

0.45 

56. 2 

7.42 x 10 

Q 

D 37 

0.46 

57. 5 

7. 59 x 10 

*48 

0.39 

43. 6 

5. 75 x 10“ 9 

oc 

y 

o 

0.39 

43.6 

5. 75 x 10“ 9 




-10 

00 

W 

0.09 

4.1 

5.42 x 10 
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2. 3 Statistical Nature of the Experiment 

As stated previously, it is not feasible to measure the precipi- 
tation reflectivity directly in all system configurations that must be 
coordinated to prevent detrimental rain- scatter interference. Thus, it 
is necessary to develop a model based upon a more common data source, 
such as a meteorological parameter that is measured generally on a 
world-wide basis. The readily available data are primarily surface 
rain measurements. 

The Virginia experiment was designed to collect data for a 
sufficient period of time to obtain a reliable statistical base, and to 
"catch" some relatively rare rainfall events in the experimental network. 
A minimum base of either rainfall or radio transmission loss statistics 
is estimated to require operation of the experiment for a period of 
at least one full year. Observation of the rare rainfall events, which 
consist of high rainfall rates that occur on the order of 1 to 5 minutes 
out of a year, appears to offer a serious challenge. However, the task 
is perhaps not as difficult as it seems because of the physical 
nature of very high rainfall rates. 

A year, which would be a totally inadequate sampling period for 
catching an extreme value of 24-hour rainfall, may be a minimum 
reasonable period for catching near -extreme one -minute rainfall rates, 
as the highest one -minute rate might be expected to fall within a year. 

The reason lies in the saturation effect at very high rainfall rates; i.e., 
there is a physical limit to the rainfall rate that may exist for very 
short periods of time. For example, the highest value of Z that has 
been reliably measured with calibrated weather radars is about 
10 mm /mm [Donaldson, 19 61 ] ; this presumably represents a near- 
instantaneous value. Assuming that this measurement is due to rain 
(more likely it was due to hail [Rinehart et al. , 19 68 3), the implied 
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rainfall rate is on the order of 1000 mm/hr. The world record one- 
minute rainfall rate was for some time thought to be 1. 23 inches (or 
1800 mm/hr) recorded at Unionville, Maryland, in 1956 [Critchfield, 

I9 60], but this record was apparently caused by sticking of the rain 
gauge in use there [ Private comm. R. K. Crane, 1970]. Credible one- 
minute records are on the order of 0. 65 to 0. 69 in. , and correspond to 
rainfall rates of roughly 1000 mm/hr. These records represent the peak rai] 
fall rates culled from hundreds of operational rain gauges throughout the 
world over a period of time of 50 years or more (one such record 
dates back to 1911). 

The rainfall rage statistical model of Rice and Holmberg [1972] 
indicates that for the Virginia area, or for Norfolk in particular, 
the number of hours of rainy minutes, T(R), at rates exceeding any 
sufficiently high rate, R, is 

T(R) = T exp (-0. 03R) hours, (18) 

where 60 T(0) = 437 is the number of rainy minutes per year of Mode 1 
(heavy, convective) rain. From (18), the rain-rate for the rainiest minute 
per average year is 204 mm/hr; for a 30-year period the rainiest 
minute should be 318 mm/hr, whereas the return period for a rate of 
1000 mm/hr is on the order of 2 x 10^ years at a median location. 

On the other hand, during the period from 24 September 
through 31 January, the 10, 000 -ft subpoint rain gauge in the Virginia 
experiment registered 1. 5 minutes of data with a rain-rate on the 
order of 100 mm/hr, and 7 seconds of rain at the rate of about 120 
mm/hr. Based on experience to this date, it seemed entirely reasonable 
to assume that, during the balance of the first full year of operation, 
one -minute rain-rates on the order of 200 mm/hr would be measured. 

A one -minute rate of approximately 150 mm/hr was measured at the 
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Poquoson tower site during August 1970, before the experiment was 
operational. The effect of an increase in rain-rate on the received 
power via rain scatter can be estimated by using an equation of the 
for m of (1 1) as 

Z = aR b . (19) 

The coefficient a is relatively unimportant in estimating the effect of 
increases in R, but the value of the exponent b is crucial. For example, 
assuming a typical value of b to be 1. 5, the effect of changing R from 
200 to 1000 mm/hr will be to increase Z, and thus the received power, 
by about 11 dB. 

The accuracy with which we can estimate the effects of the extreme 
values of R on radio systems depends on how accurately we can extra- 
polate from R « 200 mm/hr to RsslOOO mm/hr (or less). Again, assum- 
ing a in (19) to be relatively stable as a function of R and assigning the 
error in the extrapolation process to errors in the assumed value of b 
at very high R, an absolute error of 0. 25 (17%) in b, (i. e. , b = 1. 7 5 or 
1. 25 when 1. 5 was actually used) results in an error of only 1.7 5 dB in 
the calculated value of Z or transmission loss at R^J 1000 mm/hr. More 
realistically, extrapolating from R = 150 mm/hr (assuming 200 mm/hr 
is not recorded) to R = 330 mm/hr (a 4. 8-dB increase in received power) 
to construct a 30 -year distribution would involve an error under these 
assumptions of only 0.8 dB. The longer the Virginia experiment is 
operative the more likely is the probability of capturing the rare events 
(high rain-rates) and hence reduci-ng the extrapolation errors. 

2.4 Path Loss Calculations 

It is instructive to consider the fixed-beam experiment as an 
actual example of a satellite -terre strial communication systems 
problem, in which the two services share a common frequency band. 
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When frequency assignments are made on a coequal basis to these 
services, site locations must be coordinated in such a way that no 
intolerable mutual interference will result. The site separations are 
then determined on the basis of a "coordination distance" which will 
meet the established criteria. 

Procedures for calculating coordination distances are estab- 
lished internationally by Recommendation 1A, Final Acts of the Extra- 
ordinary Conference to Allocate Frequency Bands for Space Communi- 
cation Purposes (EARC) [ EARC, 1963]. Report 382 of the CCIR [ CCIR, 

1 966a] and various documents of the Xllth CCIR Plenary, New Delhi, [CCIR, 
1970] recommend modifications to these procedures, and FCC Regu- 
lation 25 [ FCC, 1966] establishes somewhat modified procedures for 
the U.S. The EARC Recommendation 1A establishes -145 dBWas the 
maximum permissible unwanted signal level at the input to an earth 
station receiver from any single terrestrial station (TS) transmitter 
with an output equal to P dBW and an antenna gain (minus feeder loss) 
equal to G dB relative to an isotropic antenna. Coordination dis- 
tances are calculated in terms of these quantities, the maximum earth 

station (ES) antenna gain (minus feeder loss) G at each azimuth and 

ES 

an earth station site shielding factor F , where applicable. This recom- 
mendation also establishes the minimum permissible sum of basic 

transmission loss and site shielding factor, L, + F , as 

b s 

H> +F s = P TS + G TS + G ES + 145 dB ' (20) 

for 4 GHz. 

Allowable site shielding factors, F , are listed in Report 382 [CCIR, 

s 

19 66a] as a function of the minimum angle of elevation of a horizon 

obstacle as seen from the earth station for each applicable azimuth. 

In current radio regulations, the maximum value of F is now 15 dB. 

s 
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In CCIR Report 382 [CCIR, 1966a], it is proposed that the maximum value 
of F should be 25 dB*. A maximum allowable value, P^g+ = 42 dBW 
in the overlapping portions of cochannel bandwidth assignments between 
a terrestrial station and an earth station,is proposed by the CCIR. 

Coordination distances at 4 GHz to the simulated earth station 
at Langley in the directions of Quantico, Norfolk, Fort Lee, and East- 
ville, respectively, are based on the relative gain of the earth station 
in the direction of these terrestrial sites, ana assume an effective 
isotropic radiated power (EIRP) of 42 dBW at a hypothetical terrestrial 
station. Using the methods outlined in Recommendation 1A of the 
EARC, the coordination distances to the Langley earth station are 
shown in table 3. 

It is necessary to add a correction factor to L^ when coordin- 
ation distances at 8 GHz are considered. These correction factors 
appear in the Recommendation 1A of the EARC [19 63] and were used 
to calculate the coordination distances at 8 GHz, which are also given 
in table 3. 

All values used in calculating the coordination distances were 
derived from the propagation curves for Zone A as given in the Final 
Acts of EARC mentioned above. The curves apply to both 4 and 8 GHz 
using the correction factors given for the higher frequency. The 
receiving antenna gain was derived from a curve especially fitted to the 
actual antenna pattern measurement [Rice et al. , 1970 ]. 

It will be noted that all transmitter-receiver combinations in the 
Virginia experiment are within the calculated coordination distance; 
therefore, it would be necessary to consider the probability of interference 

*Note: Recommendations from the International Working Group Meeting 
at Geneva, February 1971, indicate values up to 50 dB may be achieved. 
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before frequency assignments for such stations could be considered. 
A comparison between the calculated coordination distances and the 
actual great circle (GC) path distances is shown in table 3. 

Table 3. Comparison of Coordination and Actual Distances 

Transmitter Great Circle Calculated 

Site Distance from Coordination 

(TS) ES Distance 

(km) (km) 




4 GHz 

8 GHz 

Quantico (TS-1) 

178.86 

200 

1 60 

Norfolk (TS-2) 

1 6.73 

117 

100 

Ft. Lee (TS-3) 

85. 20 

120 

100 

Eastville (TS-4) 

50. 67 

125 

100 


For each of the transmitter-receiver combinations, it is also 
desirable to know the expected path loss characteristics which would 
produce a side -lobe -coupled unwanted signal. Each antenna at the 
transmitter sites will have a different gain in the direction of the 
.receiver site and the receiving antenna will have different gains in 
each of the four directions to the transmitter sites. The specific 
antenna gains were obtained by the curve -fitting procedure outlined 
by Rice et al. [1970 3. The relative antenna gains were obtained by 
considering the three-dimensional antenna pointing and great-circle 
path configurations. The difference between the vertical and hori- 
zontal polarization gains of the antennas was also considered. A computer 
program was used to calculate the relative gain factors along each 
great circle path, employing the derived antenna curves mentioned above. 
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The expected path loss for the different sites under the climatological 
conditions of Virginia was calculated using the ITS tropospheric 
prediction methods [Longley and Rice, 1968], A prediction for 
each of the great circle paths was derived, based on these methods. 

Using the above basic transmission loss models, antenna 
gain patterns, a power input to the antenna of 10 watt% and an estimated 
feeder line loss of 3 dB, the expected signal strength at the receiver 
for each of the antennas at the four transmitter sites was calculated 
using the relation 


ES 


= P 


TS 


+ g ts + G 


ES 


V 


(21) 


where L, is the path loss exceeded 99. 9% of the time. The results 
b 

are tabulated in table 4 for conditions of no site shielding. 

The tabulated signal levels are expected to be exceeded not more 
than 0. 1 percent of the time, with a probability of 0. 5. Since the receiver 
noise threshold was -130 dBm for the S-band and -125 dBm for the X-band 
channels, it can be seen that the expected signal for 0. 1 percent of the 
time thus exceeds the threshold level on all of the great-circle paths 
where the full gain of the transmitter antenna is realized. Also the off- 
path and 10, 000 ft paths at both frequencies at Norfolk are seen to pro- 
vide side-lobe interference via tropospheric propagation. The Quantico 
antenna pointed at the 20, 000 ft intersection may also provide an inter- 
ference signal about 0. 1 percent of the time. 

Similarly one can calculate the transmission loss for which 
the expected signal will be equal to the receiver sensitivity from the 
relation 


L. (Min. Sig. ) = P - Receiver Threshold + G + G 

■u To w 


TR 


ES' 


( 22 ) 
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Table 4. Calculated Values of Path Loss 
Side -lobe to Main-lobe Coupling- -Tropo. Scatter 


T ransmitter 
Site 

■ 

Path 

Common 
V olume 
(ft) 

Received 
Power 
(0. 1%) 
(dBm) 

% Time 
Rec. Sig. 
Exceeds 
Threshold 

Quantico 

f, 

A„ 

GC 

mam 

6. 0 

(TS-1) 

1 

D 

20 k 


0. 11 



E 11 ' 

30 k 


< 0.01 



F l! 

Off Path 

MSeSm 

< 0.01 


£ 

A 1 , 1 

GC 

-120.8 

0. 5 


6 

dJ 

20 k 

-144. 5 

< 0.01 



E 16 

F 16 

30 k 

-145.2 

<0.01 

Norfolk 

f. 


10 k 

H 

32. 0 

(TS-2) 

z 

F 22 

22 

Off Path 


97.0 

Ft. Lee 

f , 


GC 

1 

1.4 

(TS-3) 

3 

c 33 

10 k 


< 0.01 



D 33 

20 k 

s c hi h 

< 0.01 


f 

A 33 

GC 

-124. 3 

0. 15 


7 

c 37 

10 k 

-155. 3 

< 0.01 



D 37 

37 

20 k 

-164.5 

< 0.01 

Eastville 

f. 

A,, 

GC 

-106.9 

10.0 

(TS-4) 

4 

,44 

°44 

20 k 

-154. 3 . 

< 0.01 


f. 

B 44 

5 k 

-138.8 

< 0.01 


5 

cj! 

10 k 

-146.9 

< 0.01 




30 k 

-157.8 

< 0.01 


f 

A 45 

GC 

-121.4 

0.3 


8 

B 48 

5 k 

-152. 1 

< 0.01 



c 48 

10 k 

-160.2 

< 0.01 



E 48 

30 k 

-171. 1 

< 0.01 



48 





The percentage of time that the signal should exceed the receiver 
threshold for the calculated path loss can be read from the trans- 
mission probability curve for the respective paths. This percentage 
of time is given in the final column of table 4. These estimates show 
that the Eastville and Quantico great-circle paths at the lower frequency 
and the Norfolk paths should experience signals above threshold for 
significant periods during a typical year. These calculations utilize 
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the antenna patterns as smoothed by the curve -fitting technique; 
consequently, the actual results may show some significant deviations 
from these calculated results, based on the actual sidelobe patterns 
that existed. Also, the effect of terrain shielding has not been accounted 
for in these calculations. Data obtained from the great-circle measure- 
ments noted previously can be compared with the predictions in table 4. 

3. EXPERIMENT CONFIGURATION 

The geographic locations of the simulated* earth station (ES) 
and terrestrial stations (TS) are shown on the pictorial map of figure 1. 

A 30 -ft receiving antenna was used at the simulated ES and was aligned 
along a great circle path toward Quantico (TS-1) at a fixed elevation, 
angle of 13° 30'. Transmitting antennas at the four simulated 
terrestrial microwave radio relay stations were aligned to provide four 
common volume intersection points, at elevations along the receiver 
main beam of 5, 000, 10, 000, 20, 000, and 30, 000 ft. At Norfolk and 

Quantico an antenna was pointed "off-path" to provide the effect of a 
transmitter antenna main beam not intersecting but passing near the 
receiver antenna main beam. At the three forward stations, Quantico, 

Ft. Lee, and Eastville, a transmitter antenna was aligned along a great 
circle path toward the ES to provide a reference tropospheric propagation 
path loss measurement. 

The experiment was designed around eight continuously recording 
receiver channels each with a nominal 2 kHz bandwidth. There were five 


*The term "simulated" here is intended to imply antenna gain, sensi- 
tivity, etc., and not geometrical alignment or configuration of a 
typical satellite -terrestrial sharing system. 


27 




28 


Figure 1. Configuration of the RIPP fixed-beam experiment. 



lower frequency channels separated by 72 kHz, nominally at 3. 672 GHz. 
The three high frequency channels were nominally 7.834 GHz, 
separated by 153. 6 kHz. 

Sixteen transmitting antenna beams (nine at S-band frequencies 
and seven at the X-band) located at the four simulated TS were aligned to 
intersect the ES main beam at the common volume altitudes listed above. 
The receiver antenna was equipped with a wideband feed system 
designed to receive both the S- and X-band signals, and was aligned 
with its main beam centered in the great circle plane passing through 
Langley andQuantico. Table 5 is a complete listing of all fixed-beam 
facilities showing location, site designation, transmitter frequencies, 
antenna code, antenna purpose, and antenna pointing angle. Table 6 
lists the common volumes as computed by the conical method noted 
previously and presents all of the range parameters to the common 
volume locations, the scattering angle involved at the beam inter- 
sections, the antenna gains relative to isotropic, and the nominal 
line losses associated with each transmitter antenna. Line losses 
at the receiver are not tabulated as they were calibrated out of the 
system using the calibration technique described later. The listed 
parameters include all of those necessary to compute the K factors 
in table 1. The path designations used throughout this report 
are those of table 1. 
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Table 5. Fixed Beam Antenna Orientations 


Location 

Site 

Frequency, Hz 

Antenna 

Purpose 

Azimuth 

Elevation 

Quantico 

TS-1 

V 

3,671, 928, 000 

A n 

C. C. to ES 

152° 

13.4’ 

0° 18' 





E n 

30 k ft. 

152° 

13.4' 

3° 15.8' 





F n 

Off path 

148° 

13.4' 

0° 18’ 





D n 

20 k ft. 

15 2° 

13.4' 

1° 46. 9' 



V 

7,833, 600, 000 

A 16 

C.C. to ES 

152° 

13.4' 

0° 18' 





E 16 

30 k 

152° 

13.4' 

3° 15. 8' 





°16 

20 k 

152° 

13.4' 

1° 46. 9’ 

Norfolk 

TS-2 

f . 
2 

3,672, 000, 000 

C 22 

10 k 

337° 

17.8' 

5° 47. 5' 





F 22 

Off path 

344° 

43.4' 

0° 48' 

Ft.. Lee 

TS- 3 

V 

3,672, 072, 000 

A 33 

G.C. to ES 

101° 

6.6’ 

0° 48' 





C 33 

10 k 

93° 

43.2' 

1° 58. 9’ 





°33 

20 k 

84° 

58. 5’ 

4° 35. 9’ 



V 

7,833, 753,600 

A 37 

G. e. to ES 

101° 

8. 3' 

0° 48' 





C 37 

10 k 

93° 

43.2’ 

1° 58. 9' 





°37 

20 k 

84° 

58. 5' 

4° 35. 9’ 

Eastville 

T S-4 

V 

3,671, 856, 000 

A 44 

G.C. to ES 

237° 

58. 5' 

0° 48 ' 





°44 

20 k 

26 5° 

50.2’ 

6° 10.8’ 



V 

3,672, 144, 000 

E 45 

30 k 

276° 

59. O' 

8° 21. 5' 





C 45 

10 k 

252° 

28.1’ 

3° 14.8’ 





B 45 

5 k 

24 5° 

14.3’ 

1° 34. 1 1 



V 

7,833, 446,400 

A 48 

G.C. to ES 

237° 

58. 5' 

03 

, O 
° 





E 48 

30 k 

276° 

59. O' 

8° 21. 5' 





C 48 

10 k 

252° 

28. 1 ' 

UJ 

0 

00 





B 48 

5 k 

245° 

14.3' 

1° 34. 1' 

Langley 

ES 



Fixed 

Beam 


332° 

47.9' 

13° 14.8' 


30 
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elevation of common volume. Other symbols are defined in section 



The transmissions from each transmitter site were received at 
the ES through separate receiver channels, and thus, clear identification 
of the source of each signal can be achieved. The low and high frequency- 
antennas at each transmitter site were fixed in position, as indicated 
earlier, to point to specific locations. The transmission frequencies 
were switched at one -minute intervals in an established and fixed 
sequence, and were monitored and synchronized from the receiver 
site by means of a dedicated telephone line telemetry system. The 
transmitted power at all frequencies was held constant at a level of 
40 dBm (±0. 2 dB) by means of an electronic control system at each 
transmitter. The receiver channels were calibrated automatically 
at two -hour intervals by introducing a known signal level through 
a precision attenuator and a directional coupler mounted directly 
behind the antenna feed horn. This calibration covered a 60 -dB 
dynamic range in steps of 10 dB. This wide dynamic range was 
achieved through the use of logarithmic compression circuits in the 
IF stages of the receiver. 

The output of the eight receiver channels was recorded on both 
analog magnetic tape and a strip chart recorder. In addition to the 
received signal levels, the recorders were encoded with a time signal 
which permitted the referencing of all data to local time as well as the 
decoding of the one -minute transmitter sequence intervals. The 
strip chart was used for visual monitoring of the experiment in real 
time, for annotating the data in real time, and for back-up to the 
magnetic tape as an archival data record. 
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The one -minute antenna switching sequence was performed in a 
four -minute timing cycle, which permitted the monitoring of a total 
of 32 path combinations with the eight frequencies involved. Four 
of these combinations represent a minute during which a transmitter 
was switched to a dummy load and only the receiver noise level was 
recorded. This provided a means of monitoring the long term receiver 
performance during the experiment. Two off-path signals, one each 
from Quantico (F ) and Norfolk (F ), and six great circle paths 

X -L uu 

(A,,, A, A.., A__, A. ., and A..) were monitored. The remainder 
11 1 6 33 37 44 48 

of the available time -frequency combinations were used to monitor and 

record the sixteen paths derived from TS common volume intersections 

with the ES main beam. Some of these 32 combinations (C__, F„, A 

22 22 44 

D ) were sampled for two minutes in every four -minute cycle. The 
complete matrix of the time -frequency switching system is shown in 
table 7. 

An overall block-diagram of the fixed-beam experiment is 
shown in figure 2, in which only one of the four terrestrial transmitting 
sites is indicated. These sites were all configured in a similar 
fashion, the basic difference being the number of transmitters located 
at any one site. The latter may be seen from table 5. 

Signals that were directed over the great-circle paths in the 
configuration noted above were received in the side-lobes of the ES 
receiving antenna whenever propagation conditions permitted this 
form of interference. In order to provide a reference for this 
interference mechanism, a second great-circle path receiver 
system was installed. Three horn type antennas were located 
at the ES receiving site, directed toward the three forward 
TS sites at Fort Lee, Quantico, and Eastville. 


and 
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Figure 2. Block diagram of the RIPP fixed-beam experiment. 

























Table 7. Antenna/Path Switching Sequence 


Transmitter 

Frequency 

Antenna/Path Sequence 

Location 

^^Time 



& Site No. 





4 

Quantico TS-1 

5 

A 

11 

E n 

F 11 

D n 


f 6 

A 16 

E 16 

off 

°16 

Norfolk TS-2 

f 2 

°22 

F 22 

C 22 

E 22 

Fort Lee TS-3 

f 3 

A 33 

off 

°33 

D 33 


f 7 

A 37 

off 

C 37 

°37 

Eastville TS-4 

£ 4 

A 44 

A 44 

°44 

°44 


£ 5 

off 

E 45 

C 45 

B 45 


f 8 

A 48 

E 48 

C 48 

B 48 


Data were recorded over the same four -minute cycle (one minute for 
each path and one minute off) on an independent strip chart recorder. 
This system is also noted in the diagram of figure 2. 

A summary of the facilities and parameters associated with 
the transmit and ES receiver sites are presented, respectively, in 
tables 8 and 9. Characteristics of the great-circle receiving system 
are shown in table 9. Specific data concerning the transmitting and 
receiving systems are presented in tables 10 and 11. 

The experiment also included a rather comprehensive system 
of meteorological instrumentation indicated in a separate block of 
figure 2. Instruments to measure pressure, temperature, and relative 
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Table 8. Transmitting Site Summary 


TRANSMITTER SYSTEM 
Antenna Gains 
Antenna Pointing Accuracy 
Transmitter Output Power 
Output Power Stability 
Frequency (nominal) 
Frequency Stability 
TRANSMIT SITE 

A) Eastville (TS-4) 

B) Ft. Lee (TS-3) 

C) Quantico (TS-1) 

D) Norfolk (TS-2) 


X-BAND 

S-BAND 

41. 0 dB 

38. 8 dB 

± 15* 

± 15* 

10 Watts 

10 Watts 

+ 0. 2 dB / wk 

+ 0. 2 dB / wk 

7. 834 GHz 

3. 672 GHz 

, “10 . 

5 x 10 / day 

5x10 *0/ day 

NUMBER OF ANTENNAS 

6-FOOT DIA. 

10 -FOOT DIA. 

4 

5 

3 

3 

3 

4 

0 

2 
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Table 9- ES Receiving Site Summary 


CW SYSTEM 

Antenna Gain 
Beamwidth 

Receiver Noise Figure 
Receiver Gain Stability 
F requency Stability 

GREAT CIRCLE SYSTEM 

Antenna Gain 
Receiver Noise Figure 
Receiver Gain Stability 
F requency Stability 


X-BAND 

50.8 dB 
0.4° 

7. 5dB 
± 1 dB 

5 x 10 1^/ day 


S-BAND 

47. 5 dB 

0 . 6 ° 

8 dB 
± 1 dB 

5 x 10 _1 ° / day 


19.2 dB 
7. 5 dB 

+ 1 dB 

-10 

5X10 per day 


RADAR SYSTEM (VERLORT) 

Antenna Gain 
F requency 
Peak Output Power 
Range 

TELEMETRY 

4 leased land lines, voice grade 


Nominal 35 dB 
28 50 MHz 
200 kW 
125 nm 
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Table 10. CW Transmitting System Summary 



Frequency 

Antenna 

3. 67 2 GHz 

7. 834 GHz 

Diameter 

10 ft 

6 ft 

Gain (above isotropic) 

38. 8 dBi 

41. 0 dBi 

Beamwidth 

o 

oo 

1.5° 

Sidelobes 

-22 dB 

-21 dB 

Transmission Line 

50 ft 3. 4 dB 

50 ft 3. 1 dB 

Loss 

100 ft 6. 1 dB 

100 ft 4. 2 dB 

TWT Amplifier 



Power Output 

10 watts 

10 watts 

Gain (minimum) 

35 dB 

35 dB 

Microwave Source 


- 

Maximum Power Output 

+ 16 dBm 

+ 10 dBm 

Frequency Stability 

- 

o x 10 per day 

r ,„-10 

5x10 per day 

RF Switches 



YSWR 

< 1. 2 

< 1. 2 

Power Rating 

10 watts 

10 watts 

Life 

6 

>10 switches 

6 

>10 switches 

Power Meter 



Accuracy 

±0. 05 dB 

± 0. 05 dBm 

Stability (on 5 dBm 

± 0. 2 dB 

± 0. 2 dB 

scale used) 
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Table U. CW Receiving System Summary 



Frequency 

Antenna - 30 ft. Parab. 

3. 672 GHz 

7. 834 GHz 

Gain 

47. 5 dB 

50. 8 dB 

Beamwidth 

0. 6° 

0.4° 

Sidelobes 

-21 dB 

-18 dB 

Directional Coupler Loss 

0. 2 dB (max) 

0. 2 dB (max) 

Diplexer Loss 

0. 5 dB 

0. 4 dB 

Image Noise 

3. 0 dB 

3. 0 dB 

Tunnel Diode Amplifier 



Noise Figure 

4.9 dB 

4. 4 dB 

Gain 

20 dB 

20 dB 

Bandwidth 

~400 MHz 

~400 MHz 

Local Oscillator 

Frequency stability 

_ 1/x -10 J I 

5 x 10 per day 

-10 

5 x 10 per day 

Injection Power 

+ 7 dBm 

+ 7 dBm 

Mixer-Amplifier 



Noise Figure 

7. 5 dB 

8. 0 dB 

Gain (RF-IF) 

32 dB 

32 dB 

Bandwidth 

3 MHz 

3 MHz 

IF Center Frequency 

10 MHz 

10 MHz 

IF Amplifier 



Frequency Response 

0. 25 to 1 50 MHz 


Gain 

20 dB 


Output (Full Gain) 

1 V p-p (50 ohms) 
compression 

at 1 dB gain 

Filter, Crystal 



Pass band ripple 

t 1.5 dB 


Isolation 

70 dB 

! 


Bandwidth (3 dB) 

2 kHz 


Logarithmic IF Amplifiers 



Center Frequency 

IQ MHz 


Bandwidth 

3 MHz 


Dynamic Range 

80 dB 


Log Accuracy 

1 1 dB over 60 dB 



t 2 dB over 80 dB 


Outputs 

Video, d-c coupled 


IF, limited 
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humidity were placed at the three forward transmitting sites and 
at a meteorological site (Poquoson tower) near the ES. These meteoro- 
graph instruments recorded continuously, and were checked, cali- 
brated, and had their recording charts replaced at weekly intervals. 
Tipping bucket rain gauge instruments were located under each of the 
common volume intersection points and at the ES meteorological 
site. These instruments also recorded continuously and were 
calibrated to indicate the accumulation of each 0. 01 -inch of rain at 
rain-rates up to 7 inches per hour (177.8 mrh/hr). The rain gauges 
were capable of operating up to 45 days on a single chart record. 

This meteorological data collection system was augmented by 
radiosonde data taken twice daily at four locations surrounding the 
test site; namely. Wallops Island, Dulles Airport (near Washington, 
D.C.), Cape Hatteras, and Greensboro, North Carolina. Monthly 
weather summaries were also collected from Weather Bureau stations 
in the area. 

A primary augmentation of the above surface meteorology data 
and the radiosonde measurements was provided by an S-band surveil- 
lance radar located at the ES receiving site. This Very Long Range 
Tracker (VERLORT) radar was dedicated full time to the experiment 
and was equipped with both PPI and RHI displays which were photo- 
graphed to provide sequential records of the progress of a storm. 

A complete summary of all the subsystems in the meteorological 
configuration is presented in table 12. These data were collected 
for synoptic weather descriptions and for detailed application in 
properly identifying propagation conditions and particular interference 
mechanisms. 
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Table 12. Meteorological Subsystems 



ES*§ 

Location 

Data 

Accuracy 

Calibration 

Maintenance 

Rain 

Gauges 

5 

Poquoson 
and. common 
volume 
sub-points 

Rain 

Rate 

± 1/2% to 
20% 

Bi- 

Monthly 

Weekly 

Meteorograph 

1 

Poquoson, 
Ft. Lee, 

Eastville fr 
Quantico 

Temp. 

- 1°F 

Weekly 

Weekly 


HPigii! 

Pressure 

= C.0l' ! 

Surveillance 
Radar 
( Ve rlort) 

1 

Lang!ey 

Quai. Refi. 
(Rain) on 
PPI ii RHI 

Qualitative 

As Req'd 

As Heq'd 

Nat'l WX 

Service 

Data 

Fax Wx Maps, Radar Maps, Radar Photos Fax & Film, Rawinsonde Data 
by TTY & Fax Recorder Charts, Local WX Data Inc!. Rain Gauge Data. 


4. DATA PROCESSING AND REDUCTION PROCEDURES 
An overall data flow diagram for the Virginia fixed-beam experi- 
ment is shown in figure 3, which serves to illustrate the interface 
between meteorological and radio data as well as the detailed interfer- 
ence analysis procedure outlined in this section. All of the data records 
were collected at the Langley experiment sites and returned to the 
Institute for Telecommunication Sciences in Boulder, Colorado, for 
processing. Data were recorded in several different ways as indicated in 
figure 3. The basic radio interference data were recorded on analog 
magnetic tape, and the rain-rate (R) record was a paper chart that was 
driven by a mechanical clock on each rain gauge instrument. 

A great deal of preliminary editing was applied to these data before 
further processing took place. For example, a strip chart recording 
of all the radio signals was made at the Langley ES in parallel with the 
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magnetic tape recording. These charts as well as the rain gauge records 
were first visually scanned to determine periods of significant events. 

The latter may be periods in which radio signals were received either 
in the presence or absence of precipitation, measured precipitation 
within the configuration of the experiment, or any other significant 
weather or meteorological observation that could have an effect on the 
propagation mechanisms of interest. The times of occurrence and 
duration of these events were logged, and a listing of the data collected 
in both the radio and the meteorological systems was made. Each event 
period was compared with daily station-log reports and weekly system 
summary reports to determine the validity and accuracy of the data. 

For example, it was necessary to determine if all of the radio channels 
were in operation during the event, and if not, the outage time noted 
for each frequency. Also, each event period was given an initial class- 
ification such as "heavy convective rain, " or as a period of "non- 
precipitation. " In the latter case, the meteorological data were scanned 
for evidence of the presence of atmospheric layers or other phenomena 
that could explain any radio signals arriving via ducting and/or other 
non-precipitation scatter mechanisms. 

Time correlations were sought in the data editing process 
process between the radio, meteorological, and radar data. These 
observations served to separate convective storm periods from other 
periods (such as cyclonic rains) or to simply isolate periods of 
reception when no precipitation was observed. 
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4. 1 Radio Data 

The bistatic radio data and great-circle path data were recorded 
on a 1-inch, 14-channel, analog, magnetic tape recorder and an eight- 
channel strip chart recorder mentioned previously. The data recording 
system is outlined in figure 4, and a typical sample of a recording for 
frequency f^ is shown in figure 5. The switching sequence for the 
four-minute cycle is noted in this figure, and several significant 
precipitation interference signals are seen from the 5000 and 10, 000 ft 
common volumes during the third and fourth minutes of the timing cycle. 

A very low level signal may also be seen on the great-circle path, A . 

The significant periods of rain scatter and/or other data were 
identified and selected for digital processing, using the editing procedures 
noted above. The analog tape recorded signals for the selected periods 
were played back through an analog -to -digital (A/D) conversion 
system. All eight data tracks from the analog magnetic tapes were 
digitized simultaneously in the A/D system for the selected event periods, 
using channel multiplex and sample/hold techniques. The length of each 
digitized event period is established by the length of the longest single- 
channel data record. Thus, some channels from the analog tapes were 
digitized for periods of time when there were no significant data; a fact 
which should be kept in mind in reading the sample data tabulation 
presented later in this section. 

The pertinent factors involved in the analog playback-A/D process 
are shown in table 13. 
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RECEIVER RECEIVER 



Figure 4. Data recording diagram for the RIPP fixed-beam 
experi ment. 
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STRIP CHART RECORDER 
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Figure 5. Sample strip chart record of signals 
recorded on magnetic tape at the ES 
from the transmitter sites at Eastville 
and Fort Lee . 


46 


Signal Level - dBm 



Table 13. Factors Involved in Data Recording and Digitizing Processes 


ANALOG 

Record Speed 
Tape Playback Speed 
Speed-up Factor 


0. 6 inches per second 
15 inches per second 
25 


DIGITAL 

Record Length 
One Word 
One Character 
Two Characters 
Records per file 
Sampling Rate 


848 characters 
2 characters 
6 bits 

12 bits resolution/ word 
75 records 

16, 000 samples/second 


MULTIPLEX 

Sampling Rate _ 16000 

Speed-up factor x no. of variables 25 x 8 

= 80 channel-samples/ second 


The digitized experimental data were processed in a large scale 
digital computer (CDC 3800), and the detailed analysis of each event 
period began with the output of the computer processing. The salient 
functions and parameters that were derived from the computer programs 
were: 

1. Cumulative amplitude {in dBm) distributions for each 
clock-minute of received signal recorded at the ES on 
each frequency and for each interference transmitter 
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path (including great-circle propagation into the ES 
sidelobes ). 

2. Computation of the mean signal level in dBm, and the 
standard deviation of the clock-minute signal about the 
mean value. 

3. A count of the number of times the clock-minute signal 
crosses the median value in a positive direction during 
the minute. 

4. A tabulation of all clock-minute signal mean values, 
medians, and other deciles from the distribution of 1 
above was made from the computer analysis and recorded 
on a digital magnetic output tape for further computer 
processing outlined below. 

At the time of this report, a total of 7 0 analog magnetic tapes 
have been processed which encompasses all of the significant data 
recorded between the period of September 27, 197 0, through January 31, 

1971. The digital processing of these tapes has produced 116 archive 
data files in which the accumulated data are recorded on digital 
magnetic tape and stored. Calibration data were digitized directly 
from the analog magnetic tapes in the same manner as the radio 
signals. An example of the calibration data as stored on the digital 
archives tape is shown in table 14, which is a page of computer print- 
out for a typical calibration tabulation. The table illustrates the book- 
keeping chores performed in the data archives as each of the eight 
calibration levels for each data channel was recorded on the tape. It 
should be noted that the tabulated levels are listed in computer units 
and not units of a physical quantity. 

Following each calibration entry, the data from the distribution 
analysis were archived onto a digital tape. These records are those for each 
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Table 14. Calibration Data 
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minute of data from each channel of recording on the analog tape, and 
are cumulative values determined from 21 levels established over the 
spread of the amplitude values. A sample tabulation of these 
data for a single clock-minute is shown in table 15. In addition to the 
complete tabulation, several percentiles for each distribution are listed 
separately, as well as the mean and standard deviation values. These 
detailed descriptors of the clock-minute signals were computed and 
archived for comparison and distinction between other possible inter- 
ference mechanisms, as noted earlier in this report. All signals, includ- 
ing those recorded during periods of no precipitation, were reduced and 
analyzed in the same manner to aid in mechanism recognition and 
separation processes when required. Values in table 15 are again in 
computer units, and the listings are presented only to illustrate the 
processing, computational, and archiving details. When the appropri- 
ate calibrations are applied to these data, a tabulation of the clock- 
minute distribution percentiles (or the complete distribution) can be 
made directly in power units (dBm). An example of such a listing is 
presented in table 16, in which the data are identified for several of the 32 
possible time sequenced antenna/site/frequency combinations in the 
experiment configuration, including noise (N). The data presented 
in table 16 were taken from a period of general stratiform rain which 
occurred for approximately 5 hours on November 20, 197 0. The 
minutes which contain significant data (i. e. , mean values above minimum 
detectable levels given in section 2. 2) are indicated in the table. They 
are seen to be signals received from several great-circle paths and off- 
path beams, as well as common volume scatter data. This particular 
example serves to illustrate the compilation of data performed by the 
computer, and the read-out format from which significant clock-minutes 
of data can be recognized. 
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Table 15. Clock-minute Distribution Data 
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Table 16. Clock-minute Data Tabulation for Day 324. Levels are dBm. 
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• Denotes significant data 



In the processing of these data, a pre- and post-calibration was 
generally used with each data file. The calibration data were archived 
onto the digital tapes as noted previously. When applied to the data 
files, linear interpolation between the pre- and post- calibration was used 
whenever the particular file spanned at least half of the time between the 
two calibrations. In cases when the data file was short compared to the 
time between calibration data, the nearest calibration only was applied 
to these files. Interpolation improves the accuracy of the data by 
calibrating out system drifts and other instabilities. It is felt, however, 
that use of the nearest calibration will provide the best accuracy for 
short data files. 

The digital archives tapes were used to develop the statistical 
summaries of the tabulated clock-minute signal values as the analysis 
progressed. A second CDC-3800 program was implemented to 
compute the cumulative distributions of median signal levels for all of 
the great-circle and precipitation scatter data. The program searches 
the archived data tapes and selects all clock-minute median values that 
exceed an established threshold (such as 3 dB above noise level). These 
values were then tabulated for any or all of the 32 path combinations 
selected in program software. Any median signal level which exceeds 
the above threshold is identified with the minute at which it occurred and 
the path over which it was received. An example of a short tabulation 
of these data is given in table 17, where only and data were 
selected. These data may be selected from any single data file or 
group of files as requested in the computer program,- 
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Table 17. Clock-minute Median Tabulation 







DAY 

270 


HOUR 

ANT 

IS 

C33 

MIN 

IS 

1902 

MEDIAN 

VALUE 

-121.22 

ANT 

IS 

C37 

MIN 

IS 

1902 

MEDIAN 

VALUE 

-112.86 

ANT 

IS 

033 

MIN 

IS 

1906 

ME0IAN 

VALUE 

-123.22 

ANT 

IS 

C37 

MIN 

IS 

19 06 

MEDIAN 

V ALUE 

-117.41 

ANT 

IS 

C 33 

MIN 

IS 

1910 

MEDIAN 

VALUE 

-125.07 

ANT 

IS 

C 37 

MIN 

IS 

1910 

MEDIAN 

VALUE 

-119.66 

ANT 

IS 

C33 

MIN 

IS 

1914 

MEDIAN 

VALUE 

-123.81 

ANT 

IS 

C 37 

MIN 

IS 

1914 

MEDIAN 

VALUE 

-116.40 

ANT 

IS 

C33 

MIN 

IS 

1918 

MEDIAN 

VALUE 

-125.32 

ANT 

IS 

C37 

MIN 

IS 

1918 

MEDIAN 

VALUE 

-118.99 

ANT 

IS 

C37 

MIN 

IS 

1922 

MEDIAN 

VALUE 

-120.89 

ANT 

IS 

C33 

MIN 

IS 

1926 

MEDIAN 

VALUE 

-125.57 

ANT 

IS 

C37 

MIN 

IS 

1926 

MEDIAN 

VALUE 

-119.48 

ANT 

IS 

C37 

MIN 

IS 

1930 

MEDIAN 

VALUE 

-120.67 

ANT 

IS 

C37 

MIN 

IS 

1934 

MEDIAN 

VALUE 

-118.85 

ANT 

IS 

C 37 

MIN 

IS 

1928 

MEDIAN 

V ALUE 

-120.86 

ANT 

IS 

C37 

MIN 

IS 

1946 

MEDIAN 

VALUE 

-122.81 

ANT 

IS 

C 37 

MIN 

IS 

1958 

MEDIAN 

VALUE 

-121.62 

ANT 

IS 

C37 

MIN 

IS 

20 02 

MEDIAN 

VALUE 

-121.44 

ANT 

IS 

C 37 

MIN 

IS 

2006 

MEDIAN 

VALUE 

-122.90 

ANT 

IS 

C 37 

MIN 

IS 

2010 

MEDIAN 

VALUE 

-121.18 

ANT 

IS 

C37 

MIN 

IS 

2014 

MEDIAN 

VALUE 

-122.95 

ANT 

IS 

C37 

MIN 

IS 

2026 

MEDIAN 

VALUE 

-122.81 

ANT 

IS 

C37 

MIN 

IS 

2G30 

MEDIAN 

VALUE 

-119.37 

ANT 

IS 

C 33 

MIN 

IS 

2034 

MEDIAN 

VALUE 

-124.27 

ANT 

IS 

C37 

MIN 

IS 

2034 

MEDIAN 

V'ALUE 

-115.45 

ANT 

IS 

C 33 

MIN 

IS 

2038 

median 

VALUE 

-123.62 

ANT 

IS 

C37 

MIN 

IS 

2038 

MEDIAN 

VALUE 

-114.39 

ANT 

IS 

C 33 

MIN 

IS 

2042 

MEDIAN 

VALUE 

-123.36 

ANT 

IS 

C37 

MIN 

IS 

2042 

MEDIAN 

VALUE 

-114.39 


The program formats these selected clock-minute values into 
cumulative distributions for each path selected, based oh the percentage 
of minutes for each file or group of files. The latter can be a particular 
storm period in which the bistatic interference data were accumulated 
total of all storm periods, or a period of non-precipitation in which 
other interference mechanisms may be prevalent. Following the listing 
of table 17, the program provides a listing of the number of clock- 
minutes of data selected for each path from each data file and a similar 
listing for all previous files processed in a given run. In this manner. 
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the cumulative distributions for individual data files, particular 
storms, and/or total distributions of all previously entered data may be 
derived. The cumulative distribution functions generated in this program 
were plotted in a computer routine that uses a cathode ray tube (CRT) 
output. The data plots were recorded from the CRT on microfilm. 

An example of these function plots is shown in figure 6, printed 
directly from a frame of film. For each function, the date, starting 
data minute, and the total number of clock-minute median signal values 
accumulated in the distribution are recorded at the top of the plot. 

For long-term accumulative data, the computer also keeps a record of the 
last minute of data compiled in the function plots, and all of these 
bookkeeping notations are tabulated in the computer printout. They were 
used in later analyses to convert the percent-time scale as in figure 6 
to total minutes or hours (see section 5). 

The procedures discussed above were the primary data processing 
functions. Other parameters such as the standard deviation and the 
number of positive going median crossings within each minute were used 
when required to help distinguish and separate non-precipitation scatter 
signals. In addition, data were processed for this purpose in a 
small time series analysis system (TD-100) noted in figure 3. 

Probability distributions, correlation functions, and special analyses 
were performed in this system through hard-wired programs. Output 
was in the form of direct analog plots or photos from an internal 
CRT display, and digital magnetic tape. The latter were processed in 
the CDC-3800 system as noted in figure 3. These procedures are not 
described in detail as they were used only occasionally to augment the 
primary processing routines. The TD-100was also used to monitor J 
overall system gain and drift performance of the fixed-beam configu- 
rations through long-term analysis of noise levels and calibration data. 
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Clock-minute Median Signal Level - dBm 


ANT 02 OAT soe START Tift 1912 TOTAL TJ»€ 1176 ANT CT DAY 2?0 START Tpt 152 TOTAL Tift t’56 

. . ANT C35 Wt 270 START Tift 152 TOTAL Tift 912 



Figure 6. Cumulative distribution functions of radio clock-minute 
median signal levels for three path configurations 
and frequencies. 
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4. 2 Rain Gauge and Meteorological Data 
Meteorological data processing consisted of two distinct phases: 
processing of the rain gauge data to obtain rain-rate data for comparison 
with the radio interference data, and processing of other meteorological 
data for a synoptic description of the area weather system. These two 
phases are discussed separately. 

The rain gauge recorder charts were read manually with an 
optical comparator (a ruled- reticle equipped magnifier) yielding a 
history of the time intervals over which amounts of rain varying from 
0. 01 inch to 0.1 inch were accumulated by the rain gauge; the resolution ob- 
tainable is inversely proportional to the rain-rate, so that at rates 
over about 50 mm/hr the record was read for every 0.1 inch accumulated. 

The manually reduced rainfall record was processed with an 
XDS-940 computer program that yields both a time history of the rain- 
fall rates as recorded and a cumulative distribution of rain-rates; 
the distribution is listed by both cumulative minutes and percent of 
time for the rainfall event being reduced. An example of the 
program output is shown in table 18 in which the rain-rate (R) is 
converted to dB in the 16 log (R) column, as shown in equation (12). 

Individual rainfall events were separated into convective and 
non- convective storms (the convective class includes but is not limited 
to thunderstorms), and overall cumulative distributions were prepared for 
each type of rainfall for the entire data -taking period available. Cumu- 
lative distributions for all rainfall rates throughout the total experi- 
mental period were also prepared. In this report the total period 
examined is from September 24, 197 0, through January 31, 1971, 
inclusive. 
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Table 18. Example of Computer-Reduced Rain Gauge Data 


10, 000 ft. Subpoint, 10 Sept. 

1970, 2209- 

2229 EDT. 


Sequential Data 


Cumulative Data 

Precip. 
Rate (R) 
mm/hr. 

Time 

Int. 

min. 

16 log(R) 
dB 

Precip. 
Rate (R) 
mm/hr. 

Cumul. 

Time 

min. 

Per- 

cent 

16. 9 

0.90 

19. 7 

38. 5 

1. 98 

9. 9 

12. 7 

1. 20 

17. 7 

3 0. 8 

3. 96 

19. 9 

21. 2 

1. 44 

21. 3 

24. 2 

5. 22 

26. 2 

38. 5 

1. 98 

25. 4 

21. 2 

6. 66 

33. 4 

3 0. 8 

1. 98 

23. 8 

16. 9 

7. 56 

38. 0 

13. 4 

1. 14 

18. 1 

13. 4 

8. 70 

43. 7 

24. 2 

1. 26 

22. 2 

12. 7 

9. 90 

49. 7 

7. 5 

2. 04 

14. 0 

8. 2 

11. 76 

59. 0 

2. 5 

6. 12 

6.4 

2. 5 

19. 92 

100. 0 

Total Precip, 

. = 0.19 

inch 




Max. Rate = 

38. 5 mm/hr 





The majority of the other meteorological data, such as daily 


weather maps, radar summary maps, weather radar PPI photographs, 
and monthly weather summaries, were not processed in a routine sense. 
They were used instead as synoptic information in the separation or 
identification of mechanisms, and required analytical analysis 
using expert meteorological judgment to arrive at a conclusion. 

The radiosonde data received daily by teletype (TTY) were 
processed through an XDS-940 computer program to obtain radio 
refractive index vertical profiles, from which significant super-refractive or 
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subrefractive layers were identified and recorded in the meteorological 
log book. 

More detailed refractive index profiles were obtained for days 
during which significant radio interference signals were logged by 
obtaining copies of original radiosonde recorder charts. These charts 
were read manually and then processed through an XDS-940 computer 
program to obtain pressure, temperature, and humidity at various 
levels. This output was then processed through the same program used 
on the TTY data to obtain a more detailed radio refractive index profile. 

The meteorographs yield records of pressure, temperature, 
and humidity for one-week periods. These records were evaluated 
manually at least twice daily at the times of minimum and maximum 
temperature, and were processed with an XDS-940 computer program 
to obtain the radio refractive index at the transmitter sites. This is 
a valuable indicator of radio micro -climate. 

5. SUMMARY OF EXPERIMENTAL DATA FOR A PARTIAL YEAR 

Data from the fixed-beam experiment were collected con- 
tinuously, beginning in late September 197 0. The period of analysis 
covered in this report is from September 27, 1970, through January 31, 

1971; approximately one-third of a year. This four-month period in- 
volved slightly more than 3 000 hours of experiment time, with a total 
precipitation period of 157 hours. The latter is the time of measurable 
precipitation above 0. 1 mm/hr rates as recorded by the 10, 000-ft subpoint 
rain gauge. These data and the radio interference* data recorded during 
the precipitation periods are summarized in this section. 

*The term "interference" used here is intended to imply an unwanted signal 
from the terrestrial transmitter sites at the satellite receiver terminal. 

The degree of interference which such a signal will cause depends on many 
parameters of the satellite link, such as modulation, level of desired 
signal, etc. 
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The analysis discussed in this section is not considered to be 
conclusive, but it is useful to indicate trends that are observable 
for the period, and to illustrate the format of final results from 
the experiment. As noted in section 2, a one-year measurement period 
was considered to be the minimum required to obtain an adequate statistical 
data base from this experiment. The results presented below are the 
first cumulative summary of data that are useful to determine trends, 
and it must be empahsized again that these data represent only one-third 
of the desired year. At the time of this analysis, the experiment had 
not yet been operating during the highest rainy (thunderstorm) season of 
the year for the southern Virginia region, i.e., the summer months. 

During the precipitation period from September 27, 197 0, through 
January 31, 1971, rain scatter interference data were observed from 
the two lower common volume intersections, i. e. , 5000 ft and 10, 000 ft, 
during nearly all periods when surface precipitation was recorded. Radio 
scatter data at the common volume intersections of 20, 000 ft and 30, 000 ft 
were observed only on rare occasions from the forward trans- 
mitter sites. For example, only one clock-minute of measureable 
interference data was observed at 20, 000 ft from the cross-path 
site at Ft. Lee (path D^). This was an S-band signal that occurred 
during a 4-hour stratiform rain period in the early morning hours of 
January 23, 1971. On path D (Eastville - 20, 000 ft ) a total of 3 clock- 
minutes of measurable data were recorded, 2 minutes occurring 
during a short convective interval of a long storm on November 10-11, 

197 0. It was also during this short interval that the only interference 
data at the 30, 000 ft common volume were measured; namely 1 minute 
of data on path and 9 minutes on E . Obviously these data samples 
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are not sufficient to provide a statistically meaningful example; conse- 
quently, this report does not contain data summaries for common volumes 
above 10, 000 ft. It is anticipated that the rainy season during the summer 
months in Virginia will yield data for the common volume intersections 
above 10,000 ft £ Crane, 1972 J. 

5. 1 Comparison of Rain Scatter Data 

The experiment configuration was planned to provide redundant cover- 
age of the common volumes in both path and frequency. The redundant 
measurements provide high data reliability for the 10 , 000 and 20 , 000 ft 
common volumes, real-time data comparisons for evaluating system per- 
formance during the experiment, and a direct comparison of results 
between the S-band and X-band frequencies. 

As an example of the real-time data comparisons possible, we note 
from the switching matrix of table 7 that the 10, 000 ft common volume was 
illuminated from both the Ft. Lee and Eastville sites at both the S- and X- 
band frequ encies. during the same minute of the switching sequence. This 
permits a direct comparison of data to be made on a minute-by-minute basis 
between paths such as C and C , C and C , and between frequencies 
as vs - ^ 37 ’ Similarly, data may be compared from the 5000 ft common 
volume between the two frequencies, i.e., paths B and B Q . Scatter dia- 
gram plots of these data comparisons for several typical intervals are shown 
in figures 7 through 11. Figure 7 shows a comparison of data for paths 

and C._, taken from data intervals recorded on October 16, 1970. Note 

45 

that the diagonal line for one-to-one correspondence of these data (solid 
line) has been adjusted by 1 dB to account for the path loss difference givoi 
by the constants in table 1. The dashed line is a manual best fit to these 
data and indicates that the data samples are consistently on the order 
of 6 . 5 dB lower than data for this date. Figure 8 shows a comparison 
of these data for a storm on January 5, 1971, indicating that the data 
samples are on the order of 4 dB low for this period. 


61 



105 


o 

o 

r“H 

I 



^ap - MvaiAt axaMm-xDOTD ee o 


E 

o 

*<-i 


E 

fp 

Tt 

I 

2 

c 

w 

s 

H 

H 

D 

2 


i 

2 

O 

o 

u 

IT) 

u 


cn 

cd 

« 

T} 

c 

cd 

<u 

Cl) 


tl 

C 

cd 

1 ^ 
a) C 
Cd 

-o 

S 1 

o w 


Ih <U 

& £ 

£ -2 

go 

u > 

cd ? 
4-» O 

«} £ 

£ 

5 § 

3 o 


A! o 
8 ® 
o 2 


4) 

h 

3 

00 


62 



63 


Figure 8. Clock-minute data comparison between Ft. Lee and Eastville from 
10, 000 ft. common volume (S-band). 


However, if these two storms are combined for a minute-by-minute 
comparison of the data, the results are shown in figure 9. The regres- 
sion in this figure indicates a data discrepancy on the order of 4. 5 dB. 

The variance of these values is within the measurement accuracy of the 
experiment, which is estimated to be about ±5 dB. Individual storm 
comparisons are useful in the analysis for detecting equipment or other 
experimental problems, which will be seen in a later example. However, 
it is difficult to assess the long-range conduct of the experiment on the 
basis of such comparisons. When minute-by-minute comparisons are 

made, they cannot be conclusive because of the unknown factors such 

, * 

as storm cell size, distribution of rain within the storm, and the 
question of partial or complete filling of common volumes during each 
minute. These factors, as well as antenna orientation and system 
calibration errors^ are involved in the variations found between the 
data comparisons such as figures 7 and 8. The combination of events 
into larger statistical samples will tend to average these effects, and 
provide a better statistical norm as illustrated by the average 
regression in figure 9. 

Since the overall objective of the experiment is based upon 
statistical sampling of a long duration, it was felt that the best com- 
parison of data could be based on the total cumulative samples. Experi- 
mental error will be averaged over the longer sample, and with assurance 
that there are no events with an extreme error, the cumulative data 
will provide a basis for meaningful data comparison. Thus, the scatter 
type analysis on a minute-by-minute basis was used only to verify 
reasonable agreement between experiment measurement and rain scatter 
theory. The conclusive comparison of the data on the basis of cumula- 
tive results is presented later in this section. 

In the course of examining these data on a minute-by-minute 
basis, an interesting comparison was found for a storm which occurred 
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Figure 9. Clock-minute data comparison between Ft. Lee and Eastville from 
10, 000 ft common volume (S-band). 


on December 16, 197 0. A plot of the 10, 000 ft common volume data for 
C and C is shown in figure 10. We note ,a rather significant number 
of clock-minute means for which there is a strong signal present on 
path C^,_, while the C signal is near noise level. Regression of the 
remaining data points, however, agrees with that found in figure 9. 

The majority of the minutes during which a signal was observed for 
C and not for occurred during the morning hours, beginning around 

0900 EST and persisting until 13 00 EST, in the afternoon. A few isolated 
minutes were observed between 1400 and 17 00 EST, followed by the normal 
precipitation scatter pattern until the end of the storm at approximately 
213 0 EST. This storm represents the only event during the reporting 
period in which a serious discrepancy in the minute-by-minute compari- 
sons was observed. The storm period was analyzed for other propagation 
modes based on meteorological observations and data, and is discussed 
in some detail in Appendix C. The probability for a super-refractive 
propagation mode was found to be high on the Eastville side of the 
experiment configuration, and somewhat subrefractive on the Ft. Lee 
side. However, these conditions were not deemed adequate to permit 
the diffractive mode to explain the uncorrelated, high signal level 
minutes from the Eastville site. The transmitter antenna elevation 
angle was greater than 3° at this site. Because of the elevation angle 
and side-lobe gain considerations, superrefraction was ruled out. 
Similarly, subrefractive conditions could cause an upward bend in 

transmissions from the Ft. Lee site. Again, the elevation angle of 
o 

approximately 2 for this antenna would not produce sufficient bending 
to explain the data on this basis. 

System status for the period was carefully considered, to make 
sure the data were not due to a malfunction in a transmitter or ES 
receiver. A comparison of the minute-by-minute data at the X-band 
frequency was made (C^ vs. C^g), anc ^ * s shown in figure 11. It is 
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Figure 10, Clock-minute data comparison between Ft. Lee and Eastville from 
10, 000 ft common volume (S-band). 
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10, 000 ft common volume (X- band). 


noted that the X-band behavior was quite similar to that of the S-band 
in figure 10, with more variance in the correlated minutes. These data 
tend to confirm an anomaly in the configuration or propagation conditions 
rather than a system problem, as the S- and X-band systems were essen- 
tially independent. 

Several causes for the observed anomaly during this storm are 
possible. The combination of antenna pointing angle errors and the 
size and distribution of cells within the storm could produce the observed 
results. For example, the Ft. Lee and Eastville transmitters may not 
have been illuminating the exact same common volume. In addition, 
slight bending due to the cited refractive conditions would compound 
this problem, and minute-by-minute correlation could be destroyed. 

A second likely cause to be considered is the position of a possible 
"bright band" during the storm. If the bright band should envelop 
more of one common volume than the other (assuming that one was 
actually higher in elevation than the other), the reflectivity for the two 
paths could be significantly different. The analysis of the storm in 
Appendix C indicates that the bright band at noon (at Wallops Island) 
was slightly above 5000 ft. However, there are not sufficient data to 
confirm this altitude in the Langley vicinity. Analysis of the distribution 
and variance of the signals within the minute for the uncorrelated 
data from Eastville, indicates that these signals have the characteristics 
of rain scatter. Thus, these data were included in the total analysis, and 
it is assumed that the illuminated 10, 000-ft common volumes were not 
identical during this storm. 

As stated previously, the primary data comparisons were made 
on the basis of cumulative data for the total 4 month period analyzed 
for this report. For example, the comparison noted earlier between 
such transmis sion paths as and was made on the basis of 
the cumulative distribution data. These computer-derived distributions 
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are presented in figure 6 . Note that the separation of these distribution 
curves is on the order of 4-5 dB; a value that agrees quite well with the 
theoretical value computed from the path constants and given in table 1 

(i. e. , K - K = -124. 1 + 119. 5 = -4. 6 dB). This difference indicates 
s x 

that the clock- minute mean signal level should be 4. 6 dB higher for the 
X-band (C 37 ) than for the S-band (C 33 ) data, and is verified by- 
figure 6 . Similar comparisons were made for other path combinations 
using the cumulative distribution data presented in Appendix A. All of these 
results are summarized in the following observations: 

1) Data from path appear to be low in clock-minute mean 
value by 4. 5 dB in comparison with C 45 (which agrees also 
with figure 9). 

2) Data from paths Cgg an d C 37 agree extremely well, as 
noted above. 

3) From observations 1) and 2), we conclude that data are 
low in value by 4. 5 dB. 

4) C^g and ^ 4 g data agree well with theoretical difference of 4. 6 dB. 

5) Data from path B^g are nominally 7 dB higher than those from 

B , c . The theoretical value is 4. 5 dB from table 1. Thus, if 
4o 

B^g data are correct, the B^g data are approximately 2. 5 dB 
low in value. 

6 ) A redundant path either from Ft. Lee (TS-3) or at another 
frequency is not available for comparison directly with B^g 
data. Thus, in order to evaluate these data in a theoretical 
sense, they were compared with C^g. The two path constants 
(table 1) are the same, and it is seen from Appendix A that 
these data distributions agree very well. We conclude that 
B^g data are correct, and the low value for B^g data noted in 
5) above should be assumed. 
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data should 


7) Theoretical distributions of C ^ and 

show C generally 10. 5 dB higher than C The cumulative 
22 45 

data show C high by only 8 dB (see Appendix A), indicating 

Ct Ct 

that the data are low on the order of 2. 5 dB. 

These observations from the cumulative data for the 5000 and 

10, 000 ft common volumes show very good agreement between the 

measurements and the simplified bistatic radar scatter theory. The 

largest discrepancies are seen in the 5000 ft data from Eastville (B^) 

and the two 10, 000 ft data sets from Ft. Lee (C^ anc ^ C ). These 

discrepancies can be directly attributed to a partial blockage of 

the transmitting antennas by trees at the respective sites. This problem 

is illustrated and discussed in Appendix B. It is not clear why the 

foreground blockage appears to be 2. 5 dB higher for the S-band data 

(B ) from Eastville than the corresponding X-band (B ) data. However, 
4 d 48 

as noted in 6) above, a redundant comparison with B^g data was not 
possible, and the comparison with C^g data assumes .no height 
dependence. This may not be an absolute assumption. The 2. 5 dB 
difference in these data is, however, within the probable experimental 
error, and thus cannot be explained conclusively. 

Data from common volumes above 10, 000 ft were observed 
only sparsely during the four month period. Specific occurrences of 
contiguous samples of observed rain scatter at the higher 
altitudes were noted in the preceding section. However, isolated minutes 
of data throughout the total period were accumulated for the 20, 000 and 
30, 000 ft common volumes as well as the off-path scattered signals. 

The cumulative distribution functions for some of these data are 
included in Appendix A, even though, in most instances, the number 
of minutes in each distribution is much too small for an adequate sample. 
For instance, on the cumulative time scale in percent, a total of at 
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least 500 minutes must be accumulated before the 0. 2% probability 
level is equivalent to one clock-minute of data. Thus, the only distri- 
bution found to be somewhat significant above 10, 000 ft is that for the 

D , path, which is an "on path" common volume illuminated from the 
16 

Quantico (TS-1) site. 

The overall accuracy of the clock-minute median signal 
levels measured in the bistatic radio system has been estimated to be 
+_ 5 dB or less for all path configurations. This estimate is based 
upon a comprehensive consideration of the precision and accuracy 
of individual system components, and the physical accuracies involved 
with antenna alignments. In addition, the estimate includes those system 
effects that reflect in the use of the simplified bistatic radar equation. 
The analyses performed in arriving at this value were quite extensive 
in both aspects, and are completely outlined by Levine, et al. [19723 
for the instrumentation, and by Crane [1972] for the bistatic radar 
and raindrop distribution models. 

It should be noted that the accuracy estimate is relative to 
any one specific clock-minute of data. The dynamic error on a 
minute-by-minute basis will approach a mean-zero value within this 
estimated bound. The total experiment error however, when averaged 
over the long-term cumulative analyses presented in this section, will 
be an rms value of the minute-by-minute errors, and thus considerably 
smaller than the + 5 dB bound. 
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5. 2 Comparison of Radio and Surface Rain Rate Data 

Sections 2. 1 and 2. 2 have outlined the methods of comparing data 

in the Virginia experiment. It was noted that equations (15) and (16) 

can be used directly to compare measured values of interference with 

computed values where the surface rain rate (R ) is substituted for 

s 

the unknown rate (R^) in the common volume. Also, equation (17) can 
be used to compute an effective rain rate (R g ) from the measured 
values of interference power (P r ) over any path. Both of these 
methods have been applied to the experimental data, and the results 
are presented in this section. 

The analysis involving equations (15) and (16) begins with the 
distribution of surface rain-rates measured in the experiment. These 
distributions for the entire period from September 24, 1970, through 
January 31, 1971, are given in Appendix A. They include the distri- 
butions of the total rain measured at the 10,000, 20,000, and 30,000 ft 
subpoint gauges. The convective rain-rates measured at the 10,000 ft 
gauge are generally higher and longer in duration than those measured 
at the other points in the configuration. Appendix A also shows cumula- 
tive distribution curves (figure A-4) that combine the measurements 
from all three subpoint gauges, together with a summary table for all 
rainfall measured in the region of the experiment. The highest rate 
of 120 mm/ hr noted in these distributions was measured prior to Septem- 
ber 27, 1970, at which time the precipitation scatter measurements 
were begun. Thus, in the comparative analysis presented below, the 
rain-rate data used do not include the period prior to September 27, 1970. 
The highest rate recorded during the bistatic radio measurements, 
is on the order of 90 mm/ hr, as seen later in this section. 
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Data from the bistatic CW radio system have been analyzed 
for each clock-minute in terms of median signal level. However, the 
mean signal power should be used in the bistatic radar equation. Thus, 

before the data were plotted in the comparisons between computed and 
observed values, the median signal level was converted to a mean power 
value for each clock-minute. The conversion was made as follows: 

1) In detailed analysis, the distribution of the received signal 

within the minute (as recorded from the ES receiver) was 

* 

found to be nearly Gaussian in dBm [Hull et al. , 1970] . 

2) The normal distribution in dBm at the receiver output means 
that the power in the input signal is distributed as a log- 
normal variable. The mean power of the log-normal input 
signal can then be computed from [Beckmann, 1967], 


P = 



(23) 


where p = the mean power of the log-normal signal in mw 

P m = median power of the normal distribution in dBm 
o = standard deviation of the normal distribution in dB 
c = constant =10 log^e = 4. 3429 


Converting (23) to dB, we have 


P= 10 log p = O 2 / 2 + P m 


P m + 


1. 5 dB 


(24) 


*Hull, J. A. , P. Li. Rice, R.W. Hubbard, and G.D. Thayer [ 19703, An ex- 
periment for the study of earth station-terrestrial station interference, Pt. 1 
Interference model for rain scatter and fixed-beam experiment, Telecom- 
munications Tech. Memo. OT/lTSTM 14, unpublished report. 
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where the 1. 5 dB is estimated from the average variance observed for 
the normal distribution of rain-scatter signals within the minute. 

The ES receiver process included a linear-log detector and a low-pass 
filter at the data output. This processing of a signal log normally- 
distributed in power, results in a slight reduction in signal power at 
the receiver output. This reduction is estimated to be 1 dB for the ES 
receiving system used. Therefore, in all of the data plots presented in 
this section, each clock-minute median value has been corrected by 
+2. 5 dB in order to report clock-minute mean power levels. This same 
correction factor has been used by Crane [l972 ], but is derived in a 
more rigorous fashion, starting with the premise that bistatic radar 
data are Rayleigh distributed. 

In the analyses of these data, the measured surface rain-rate 
was converted to an effective or predicted received power level, using eithe 
equation (15) or (16) and the constants in table 1. A direct comparison 
of the converted rain-rate statistics with those measured in the bistatic 
CW radio system at S- and X-band frequencies is shown in figures 12 
through 15 for typical path intersections at the 5000 and 10,000 ft 
common volumes. In each case, the 10, 000 ft subpoint rain gauge data 
were used as noted above, deleting rain data for those periods of time 
when the radio signals were lost due to equipment or power failures. 

The total data sample (approximately 4 months) has not been sufficient 
to determine how well the bistatic radio data will relate to surface values 
at the higher rain-rates. However, significant trends may be observed 
in these data presentations, and in the following analyses. 

For most engineering purposes, the intercomparison of rain- 
rate and bistatic radio data (as seen in figures 12-15) indicates that 
an adequate prediction of the probability of precipitation scatter can be 
made from the distribution of surface rain-rates for common volumes 
up to at least 10,000 ft altitude. This is particularly true if adjustments 
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common volumes (S-band data). 





Run No. 15 (cumulative) 
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common volumes (X-band data). 
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NUMBER OF MINUTES 

Figure 14. Intercomparison of data from Ft. Lee (TS-3) for the 10,000 ft common volume 
(S-band data). 



Run No s. 1,4 through 12 
27 Sept., 1970 -31 Jon., 1971 
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Figure 15. Intercomparison of data from Ft. Lee (TS-3) for the 10, 000 ft common volume 
(X-band data). 




to the data for paths (figure 12), (figure 14), and C (figure 15) 

are made in accord with the observations in section 5. 1 for the antenna 

foreground blockage at the transmitter sites. For example, the data in 

the solid curves in these figures have not been adjusted. However, the 

noted foreground factors have been applied to the above three paths, 

and the corrected rain- scatter distributions are shown by the broken- 

line curves. The latter are seen to be very close to the predicted 

distribution based on the surface rain measurements. 

These data plots are displayed in dBm, and are therefore 

relatively insensitive at the higher rain-rates because of the logarithmic 

scale. A more sensitive data comparison can be made by using (17), in 

which an effective rain-rate, R , is calculated from the measured 

e 

interference power and compared with the surface value, R^, on a 

linear scale. To illustrate the sensitivity of these calculations, the 

following analyses \yere made. First, we wished to observe the data for a 

single storm. Figure 16 is an analysis of a typical stratiform rain 

period which occurred on November 10-11, 197 0. It can be seen that 

the radio data used to calculate an effective rain- rate, R , for this 

e 

period were at best only a gross indication of the measured surface value, 

R . In addition, the highest measured rain-rate at the surface was 
s 

not observed in the radio data, since it was measured for only 1 minute 
at the surface. There is no minute-by-minute guarantee that a surface 
rate will be observed in any given common volume aloft. A typical short 
convective storm was also analyzed which lasted for a period of about 
105 minutes. It was found (data not presented) that the radio data were 
fairly accurate in predicting the observed rain-rates below about 22 mm/hr, 
but data were lacking for comparison above this rate to the maximum 
of 90 mm/hr which lasted less than a minute. These results are 
considered typical for short duration storms, particularly where the 
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Figure 16. Comparison of measured, and calculated rain-rates for a single stratiform rainstorm. 


radio data are not collected continuously on any one path. 

The primary reason for only partial agreement in these data 
calculations is the relatively small sample size of a single storm and 
the possible resolution of the rain measurement and the bistatic radio 
data. Better agreement is expected between R and R when compari- 
sons are made over longer statistical samples, from which the radio 
data may be resolved to relatively lower intervals of time. Examples 
of such comparisons are shown in figures 17, 18, and 19, which cover 
all of the accumulated data for the period through January 31, 1971. 

The data presented in figure 17 are complete for the entire period, while 
those plotted in figures 18 and 19 are complete except for short periods 
of data during the run numbers that are omitted in the figures. The 
missing data were lost due to transmitter failures and a station power 
failure at the Eastville site (TS-4) during a storm in September 197 0. 

The rain-rate distribution has been computed accordingly in these 
figures, by deleting the storm periods corresponding to the noted 
missing data. As in the previous figures, the data from paths with an 
antenna foreground blockage have been plotted directly as measured, 
and for corrected values based on the factors given in section 5. 1. 

Figures 17, 18, and 19 demonstrate that for cumulative data 
over many storms, better agreement between the surface and effective 
rain-rates is achieved. They also indicate that the R curves will depict 
the higher rain-rates aloft more accurately over longer periods of time. 
The total data sample reported is small relative to the year-long statistics 
sought in the total experiment, but results of this four month phase 
have been found to be informative. It is perhaps possible to detect 
in these data some height dependence between the 5000 ft and 10, 000 ft 
common volumes. For example, note in both figures 17 and 18 that the 
5000 ft data are lower than the 10, 000 ft data as the curves approach 
the high rain-rates. If this trend persists in later accumulations, it could 


82 



RAIN: TOTAL (MIXED) 
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Figure 17. Comparison of measured and calculated rain-rates for total cumulative rainfall. 



RUN NO. 15 (CUMULATIVE) 
27 SEPT. ,1970- 31 JAN., I! 
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Figure 18. Comparison of measured and calculated rain-rates for total cumulative rainfall 
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NUMBER OF MINUTES 

Figure 19. Comparison of measured and calculated rain-rates for total cumulative rainfall. 




imply that the heavier rain is more dispersive in nature at 5000 ft as 
compared with rain at higher altitudes. This conjecture does not seem 
too improbable based on a study of highly convective thunderstorms 
[Culnan et al. , 1965], but current data such as those in figures 17 and 
18 are certainly inconclusive. Some empirical check on height depen- 
dence of this type, particularly above 10, 000 ft, is expected as a final 
result of the fixed-beam experiment [Crane, 1972]. 

In all of the comparative analysis presented in this section, 
the clock-minute median signal levels and the number of accumulated 
minutes have been determined from the computer-derived cumulative 
distribution functions mentioned previously. As the data were processed, 
two distribution functions were compiled and plotted on the CRT at the end 
of each run. These were the distributions for the current run and a total 
cumulative plot for this and all previous runs. Fifteen computer runs 
were made for the analysis presented here, and the complete cumulative 
plots for the data analyzed are given in Appendix A. The starting day and 
minute at which data were first assimilated, and the total number of 
minutes accumulated for the plot are printed at the top of each computer 
printout. This information was compiled by the computer for each path, 
and was written automatically on the CRT plots. The total minute count 
was used to convert the percentage scale on the abscissa to total minutes 
as seen in the data plots in this section. The computer counts do not 
reflect calibration periods, however, and some adjustment was necessary 
in the analysis. For example, the calibration period was 15 minutes (usual! 
every two hours) during which time data were not recorded. These 
minutes were accounted for manually in the analysis to determine the total 
minutes for each distribution, using the following criteria. If data were 
observed both before and after a calibration period, the 15 minutes of 
calibration time was added to the computer count assuming that data 
would have been recorded. If data were observed only before or after a 
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calibration, the calibration minutes were not added to the computer count. 
A good agreement between total number of rainy minutes and the 
bistatic radio clock-minutes has been found using this procedure. 

Evidence of this agreement is seen in the data plots in this section. 

The total number of bistatic radio clock-minutes compiled 
in the cumulative data given in Appendix A were calculated in accordance 
with the switching sequence given in table 7. Note from this table that, 
in most instances, the data were sampled for only one minute out of each 
four -minute time block. Thus, the total number of minutes of radio 
data used for comparative purposes in this section was derived by multi- 
plying the number of accumulated minutes by 4. This, of course, 
assumes that the distribution for each unobserved minute within the block 
would have been essentially the same as the observed minute. As a check 
on this assumption, the data for path C are informative. This path was 

U La 

monitored on the first and third minute of each four -minute block. The 
data for the first and third minute were accumulated separately in the 
computer routine, and, as can be seen from the cumulative functions in 
Appendix A, the two distributions are nearly identical. This provides 
reasonable assurance that the sampling rate for the radio clock-minute 
data was adequate. 
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5. 3 Relationship to Rainfall Statistics Model 

The statistical model of surface rainfall rates developed by Rice 
and Holmberg £l97 2] is based upon the sum of individual exponential 
modes of rainfall rates, each with a characteristic average rate, R. 

The temporal statistical data used in developing this model were taken 
from a great many rain gauge stations and for periods spanning many 
years of data. The data have shown that the average rainfall rate for 
each of the modes is fairly constant, but the coefficients that represent 
the total number of rainy t-minute periods for each mode are relatively 
more variable from year to year and station to station or climatic 
region. According to the descriptive analysis proposed for summarizing 
the cumulative time statistics of surface rainfall rates at a point. 

Rainfall = Mode 1 rain + Mode 2 rain, (25) 

where the exponential distribution chosen to describe "Mode 1 rain" 
corresponds to a physical analysis of thunderstorms, while "Mode 2 
irain", represented by the sum of two exponential distributions, is 
simply "all other rain". 

The average annual total rainfall depth, M, is the sum of the con- 
tributions of the two modes, M^ and M^; 

M = Mj + M 2 mm, (26) 

and the ratio of "thunderstorm rain", M^, to total rain, M, is 

P = M V M • (27) 

The number of hours of rainy t-minute periods for which a surface 
point rainfall rate R is exceeded is the sum of contributions from the 
two discrete modes: 
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T t <R) = T lt q lt (R) + T 2t q 2t (R) hours, (28) 

where 

T lt = M 1 /R lt = 0 M/R lt hours, (29) 

T 2t = M 2 /R 2t = (1-0) M/R 2t hours, (30) 

and 

q lt (R) = exp (-R/R lt ), (31) 

q 2t ( R ) = 0. 35 exp (-0. 453074 R/R 2t ) 

+ 0. 65 exp (-2. 857143 R/R 2t ), (32) 

The quantities R^. and R 2t are average annual Mode 1 and Mode 2 
rainfall rates. For t = 1, the parameter R represents clock-minute 
rates, and the values from the model are R^ = 100/ 3 mm/hr, and 
R 2 j = 1. 75505 mm/ hr. Therefore, for clock-minute rates, T^(R) 
may be written as 

Tj(R) = M jo. 03 0 exp (-0. 03 R) 

+ 0. 2 (1 -0) £exp (-0. 258 R) + 1. 86 exp ( - 1 . 63 R )j j- hours 

(33) 


Figure 20 shows T^(R)/M plotted versus R and the parameter 
j8, and Figure 21 and 22 are world maps of M and j8, respectively. 

In general, R^ t and R 2 ^ are functions not only of t, but also of 
M and the number of rainy days, D/ 24, in an average year. Figure 23 
shows average-year cumulative distributions T^ (R) versus R for 
t = 1, 5, 30, 60, 360, and 1440 minutes (1 day) for M = 1000 mm (40 inches). 
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Figure 20. Normalized cumulative time distributions of clock-minute surface rainfall rates. 
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Figure 22. World contour map of the ratio of "thunderstorm" rainfall to total rainfall depth. 
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TOTAL NUMBER OF HOURS, T t (R), OF RAINY t -MINUTE PERIODS FOR WHICH R IS EXCEEDED 
Figure 23. Average-year cumulative time distributions of clock-t-minute surface rainfall rates. 




D = 3M + 24 hours (a good average for U.S. stations), and $ = 0 . 125. 

For t = 60 minutes, Figure 24 shows similar estimates for Norfolk, Va. , 
compared with the average of observed clock-hour rate durations T^q(R) 
for the period 1951-1960. The cumulative distribution is shown for 
both R(T 6q ) and 16 log 1Q R(T 6q ). 

The data upon which this model is based include; (1) average-year 
cumulative distributions of clock-hour rates (mostly for the years 1951- 
1960), for a total of 63 stations, 49 in the continental U.S,*, (2) distribu- 
tions for all 15-year averages with t/ 60 = 6, 12, and 24 hours for 22 of 
these stations; (3) short-duration excessive precipitation storm data 
with t = 5, 10, 15, 20, 30, 45, 60, 80, 100, 120, and 180 minutes; 

(4) a U. S. map of the highest 5 -minute rates observed in a two year 
period; and (5) maximum monthly accumulations for the period 1931- 
1960 for 17 U, S, stations and 135 additional stations reported by the 
World Meteorological Organization. 

The data sample that has been analyzed for this report is not 
adequate to make a conclusive comparison with the rain-rate model. 

The model is based on conditions for an average year, and the measured 
rain data from the experiment represents approximately 35% of the year. 
However, the statistical rain-rate data presented in Appendix A have 
been used to make a comparison with the model for a partial year in the 
southern Virginia region. For example, the statistical model given by (33) 
for this region estimates a total of 7. 34 hours of rainy minutes of Mode 1, 
and 513. 63 hours for Mode 2 in an average year. Combining the two 
modes, the model estimates a total of 521 hours of rainy minutes. It is not 
strictly correct to adjust this estimate for a partial year, as the distri- 
bution of rainfall over the year is not uniform. However, if this is done, 
the model would estimate the total number of hours of rainy minutes for 
the experiment period (35% of a year) to be 182. This compares favorably 
with the actual measured value of 157 hours for the 10, 000 ft subpoint 
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Figure 24. Average-year cumulative distribution of clock-hour surface rainfall rates for Norfolk, 
Virginia. 




rain-gauge. The accuracy of the gauge method to resolve rain-rates 
of 1 mm/hr or less is, however, low. Comparing only those rates 
above 1 mm/ hr, the model estimates 76 hours of rainy-minutes for the 
partial year, and the measured value is approximately 61 hours. 

In these calculations from the model (33), the values of M = 1146mm 
and = 0. 2134 for Norfolk, Va. have been used as shown in figure 24. 

The results are seen to compare quite well, and provide the only point 
check possible for the Rice and Holmberg model at this time. 


96 



6. SUMMARY 


This report has described an experiment designed to obtain 
a statistical measure of bistatic radio interference between a simulated 
satellite receiver and terrestrial microwave communication systems. 

The measurements were made in the 4 and 8 GHz frequency bands which 
these services share in assignment. One of the basic objectives of the 
experiment is to determine how well measurements and/ or known tempo- 

£■ 

ral distributions of surface rainfall can be used to predict interference 
between these services due to precipitation scatter. Rain or other hydro- 
meteors can cause serious interference levels by scattering energy 
from one antenna system into another. Such interference is more probable 
when hydrometeors exist in the space defined by antenna beam inter- 
sections, but can also occur when energy is scattered from one antenna 
main beam into the side -lobes of another. 

Methods for prediction of the precipitation scatter interference 
are necessary to establish new coordination procedures in the planning 
phases for satellite and/ or terrestrial microwave communication faci- 
lities. Previous recommendations of the CCIR for coordinating between 
these services did not include specific consideration of precipitation 
scatter interference. The Virginia experiment and others looking into 
this problem have shown that this type of interference can be significant 
when antenna beams are likely to intersect at altitudes below about 
20,000 ft. 

Results of the experiment presented in this report cover the 
period from mid-September, 1970, through January, 1971. This repre- 
sents approximately one-third of the planned experiment year, and does 
not include the high-rain season for southern Virginia. The results for 
the period, however, have been very useful in testing the stated prediction 
hypothesis, and have been used by the CCIR and the World Administrative 
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Radio Conference (WARC) in their consideration of rain- scatter models 
for new coordination procedures. The separation of transmitter and 
receiver sites in the experiment are shown to be within or near the co- 
ordination distances previously recommended for the satellite- 
terrestrial services. Thus, the observed duration and levels of inter- 
ference due to the precipitation scatter mechanism were found to be 
quite useful in developing the proposed new coordination procedures. 

Measured sur-face rain-rates and the bistatic radio interfe- 
rence levels observed are directly compared in the data analysis. 
Temporal distributions of the surface rain-rates are converted to 
equivalent radio interference levels using the simplified bistatic radar 
theory, and a fairly standard relationship between Z (radar reflectivity) 
and R, the surface rain-rate. This conversion derives a predicted inter- 
ference distribution curve for each common volume intersection in the 
experiment, and is compared directly with the distribution of the mea- 
sured radio interference. The process is also reversed in the analysis, 
converting the bistatic interference level to an effective rain-rate 
distribution in the common volume. These curves are then compared 
with the temporal distribution of the surface rain-rate, measured by 
tipping -bucket gauges located at the sub-points of the common volume 
intersections of the antenna patterns. 

The first comparison has shown that the surface rain-rate 
calculations could have been used to predict the clock-minute interference 
signal distribution to within 1 to 3 dB for the common volumes up to 
10, 000 ft above the surface. The data sample was insufficient to test the 
prediction above this altitude. The second data comparison has shown 
that the calculated effective rain- rate distributions in the common 
volumes match the measured surface distribution to within ±5 mm/ hr 
or better for surface rates below about 40 mm/ hr. Again, the data 
sample and resolution of the bistatic radio data preclude a comparison 
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above this rate. However, the data plots indicate that the radio data will 
reasonably predict the higher rain-rate distribution when the data sample 
has been increased. 

In addition to the comparison of the measured temporal 
statistics noted above, a cursory comparison of the experiment results 
with a model of the temporal rainfall statistics developed by Rice and 
Holmberg [l972] has been made (section 5. 3). Even though the total data 
sample is less than the desired year, a reasonable agreement between 
the observation and the model has been seen. This aspect is important, 
as the Virginia experiment can be used as an "anchor point" in geographic 
area to test the Rice and Holmberg or other models of the temporal rain- 
fall distribution. Successful comparison of the model and measurement 
for the Virginia region then permits the extension of the model for use 
in predicting interference in other geographic regions. 

In summary, it is concluded that the Virginia precipitation 
scatter experiment has provided the first well-defined temporal data 
relative to the rain-scatter interference problem between satellite and 
terrestrial communication systems. This experiment, together with 
efforts to model the problem, has been instrumental in defining new 
coordination procedures which will enhance interference -free operation 
of tnese services. 

As noted in the Foreword, the balance of data from this experi- 
ment (February through September, 1971) is on file at the ITS Laboratories . 
Available funding did not permit the continued computer analyses. 

A hand analysis has been made, however, for a limited number of paths 
at the S-band frequencies, and are reported by Crane (1972). His results 
include the high-rain season (the summer months) for Virginia, and sub- 
stantiate the conclusion drawn in this report for rain- scatter statistics 
up to 10, 000 -ft altitudes. The additional data also reflect an observed 
height dependence of the scatter effects above 10, 000 ft, as noted in 
section 5.2. 
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APPENDIX A 


Total Cumulative Distribution Plots of 
Surface Rain-Rate and Bistatic Radio Data 

This appendix contains the total cumulative distribution plots 
of the surface rain-rate measurements, and the bistatic radio data 
accumulated in the computer process described in section 4. 1. A 
complete tabulation of the cumulative rain-rates used for the data com- 
parisons in section 5 is also given in table A-l. A total of 15 computer 
runs were used in analyzing the radio data. Table A-l indicates the 
rain tabulation relevant to each of the computer run numbers. 

Figures A-l through A-4 present the cumulative distributions 
of the surface rain-rates, measured at the 10,000, 20,000, and 30,000 ft 
subpoint gauges.* Figures A-5 through A-15 present the cumulative 
distribution functions of the bistatic radio data for the 5000 and 
10, 000 ft common volumes used in the comparative analysis in section 5, 
and for other path configurations where data were accumulated. 


The data for the 5000 ft subpoint gauge are not included, as this gauge was 
out of service for an extended period due to malfunction. 
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Figure A_6. Cumulative distribution of B^g data. 
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Figure A-7. 


Cumulative distribution of C^ data. 
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jure A-10. Cumulative distribution of C Ar . data. 
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Figure A-12. Cumulative distributions for the 20, 000 ft 
common volume data. 
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Figure A-13. Cumulative distribution of E^g data. 
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Figure A-14. Cumulative distributions for the off -path data, F 22 * 
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Figure A-15. Cumulative distributions of the Great Circle 
path data. 
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APPENDIX B 


Transmitter Foreground and Horizon Effects 

In section 5 . 1 it was determined that a few of the bistatic CW 
radio signals, based on clock-minute distributions, were lower in value 
than predicted by the bistatic theory. On a minute-by-minute basis 
the individual signal means will fluctuate around the predicted value due 
to variations in the actual amount of rainfall within the common 
volume. The calculated expected value is based upon the assumption 
cf complete uniform filling of the common volumes by rain. Analysis 
of the clock-minute data for all of the common volumes and paths over 
the total data sample have produced the cumulative distributions in 
Appendix A. From the comparison of these curves, it was deduced that 
the data from the B^, C and Cgy path/ frequency elements of the 
experiment configuration were low. An antenna foreground effect at the 
transmitter sites was cited as the most probable explanation. 

Figure B-l is a photograph of the foreground and horizon tree- 
line in front of the transmitting antennas located at the Fort Lee (TS- 3 ) 
site. The small crosses indicate the tops of the tree lines or the 
radio horizon for each path configuration, including the great circle 
(GC). The circles illustrate the approximate location and size of the 
transmitter beam pattern at the horizon for the S-band system. It can 
be seen that a shielding effect did exist on the 10,000 ft orientation, 
as the lower portion of the transmitted beam directed to this common 
volume was in the horizon tree-line. It was seen in section 5 . 2 that a 
correction of 4 . 5 dB brings the C33 and C37 data into line with the other 
data. 

It is noted from figure B-l that the great circle data from Ft. Lee 
was also affected by trees in the foreground. 
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Figure B-l. Foreground and horizon for Ft. Lee 
transmitter site (TS-3). 


The oth^r configurations from the site, however, did not require 
allowances for the foreground. The methods available for calculating 
the average foreground effects provide estimates consistent with those 
observed, but cannot be conclusive without supporting field measure- 
ments to provide good estimates for individual cases. 

Figure B-2 presents two photographs taken at the Eastville 
(TS-4) transmitter site. Here we note that a significant foreground 
effect was possible for the 5000 ft common volume orientation, which 

is seen in the data for paths B and B„„ as noted in section 5.1. 
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Figure B-2. Foreground and horizon of the Eastville transmitter 
site (TS-4). 
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APPENDIX C 


Example of Meteorological Analysis 

The following synopsis of a storm that occurred in the experi- 
ment area on December 16, 1970, is presented as an example of the 
general weather and meteorological analysis applied in the Virginia 
experiment. This particular storm was selected to help explain the 
anomalous data noted in section 5. 1. 

At 0700 EST, a low pressure system was centered over Cairo, 

111. , with an occluded front extending to SE Alabama. From this point, 
a cold front extended SW into the Gulf of Mexico, and a warm front 
extended ENE to the Atlantic coast near Charleston, S. C. , and SE into 
the Atlantic ocean. A shield of light to moderate rain covered most of 
the SE states, with thunderstorms along the cold front. This storm and 
its associated frontal system moved rapidly to the ENE, the low pressure 
center reaching a point over New York City by 0700 EST on the 17th. 

Widespread warm-frontal, stratiform type rain occurred with 
this system over the experimental area from about 0900 EST on the 
l6th until passage of the cold (occluded) front through the area at approxi- 
mately 2130 EST. The surface maps for these two dates are shown in 
figures C-l and C-2. 

Cloud tops associated with this frontal system were detected by 
weather radar and found to be generally between 16,000 and 24,000 ft 
over the experiment area. Rainfall began at about 0800 EST on 
December 16 in southern Virginia and about 1100 to 1200 EST in northern 
Virginia, generally ending at 2300 EST as reported by National Weather 
Service stations in the area. Rainfall was fairly uniform, averaging 
1. 24 inches at 6 stations in the experiment area and ranging from 0. 98 
inch at Dulles International Airport (Sterling, Va. ) to 1. 76 inches at 
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Norfolk, Va. Maximum hourly amounts averaged 0. 30 inch (7. 6 mm/ hr) 
and ranged from 0. 24 inch at Wallops Island, Va. and Washington 
National Airport, to 0. 39 inch at Norfolk, Va. The rain gauge 
network in the experiment reported average rainfall of 1. 57 inches, 
ranging from 1. 43 inches at the 5000 ft subpoint to 1. 79 inches at the 
10,000 ft subpoint. Maximum hourly amounts (not clock-hours as the 
National Weather Service stations report) averaged 0. 42 inch, ranging from 
0. 32 inch at the 5000 ft subpoint to 0. 50 inch at the 30, 000 ft subpoint. 

The reported rainfall and maximum hourly amounts from the experiment 
network compare favorably with the NWS report for Norfolk, Va. , 
the closest reporting station. 

The meteorological network indicated that cold frontal passage 
occurred from a WSW direction at approximately 2000 EST at Ft. Lee, 

2100 EST at Quantico, 2145 at Poquoson, and 2230 at Eastville. 

A brief heavy shower, apparently of convective origin due to 
lifting just prior to passage of the cold front, occurred at the 30, 000 ft 
subpoint at approximately 2100 EST. Duration of this shower, which 
was imbedded in stratiform rain, was 4. 8 minutes, during which time 
0. 20 inch of rain was recorded - a rate of 65 mm/ hr. The maximum 
one-minute rain-rate with this shower was 85 mm/ hr. 

The PPI display of the VERLORT radar located at the Langley 
site was photographed between 0915 EST and 2250 EST on the 16th. 

During ’this time the radar indicated widespread rainfall of varying 
intensity throughout the experiment area in patterns consistent with 
stratiform rain. Especially strong echoes were photographed between 
1900 EST and 2200 EST. The frontal passage was indicated to have 
occurred between 2134 EST and 2143 EST along a line generally parallel 
to the Langley- Quantico earth station beam, with the front moving in a 
direction perpendicular to the earth station beam. This analysis agrees 
with the rain gauge network data, which indicate virtually simultaneous 
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cessation of rain at all sites, the ending times ranging from 2135 EST 
at the 10, 000 ft subpoint to 2140 EST at the 30, 000 ft subpoint. 

Rawinsonde data for the date showed superrefractive conditions 
to be prevalent over the eastern part of the experiment area until at 
least noon, but disappearing by 1800 EST. Quick-look significant level 
data received by teletype showed only some elevated layers that would 
have had little effect on the radio system. However, detailed 
refractive index profiles derived from the rawinsonde recorder charts 
yielded different results. Surface ducts were found at Wallops Island, 
Va. , Hatteras, and Greensboro, N. C., but not at Washington, D. C., 
for 0630 EST. The temperature, relative humidity and N profiles for 
Wallops Island are shown in figures C-3 and C-4. At approximately 
1230 EST, a special rawinsonde was launched at each of those four 
stations. These flights indicated superrefractive (but not ducting) at 
Wallops Island, Hatteras, and Greensboro, with Washington being 
subrefr active. At 1830 EST rawinsonde at all four stations showed 
saturated profiles and slightly subnormal refractive index gradients. 

There is thus a good probability, based on this analysis, that 
radio signals could have been effected by a superrefractive/ diffraction 
mode during the morning hours of December 16, with the probability 
of superrefraction on the Eastville paths being high, and the probability 
for Ft. Lee and Quantico being low. There is also some evidence that 
subrefractive conditions could have been prevalent on the Ft. Lee side 
of the configuration. These observations, including the description of 
the varying intensity of the associated rainfall have been used to help 
explain the anomolous radio data observed on this date. 
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Figure C-4. N-Profile. Wallops Island, Virginia. 
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